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INVESTIGATION  ON  A  TWO-PHASE 
PROBLEM  IN  CLOSED  PIPES 

by 

Walter  Hans  Graf 

I.  ABSTRACT 

This  investigation  is  a  contribution  to  our  knowledge  of  the  two- 
phase  flow  problem .  Two  measuring  devices  for  water -sand  mixtures 
were  tested,  a  Loop  system  and  a  special  Venturi  meter.  The  Loop 
system  consisted  of  two  identical  vertical  pipe  sections  with,  opposite 
flow  direction .  Head  loss  readings  were  obtained  at  both  pipes .  The 
summation  of  the  readings  could  be  correlated  with  the. flow  rate  and 
their  difference  appeared  to  be  proportional  to  the  sediment  concentra¬ 
tion.  The  experimental  results  were  obtained  in  a  piece  of  equipment, 
especially  developed  for  the  study  of  water -solids  mixtures  in  pipes. 

In  the  second  part  a  Venturi  meter  was  investigated.  The  Venturi  meter 
was  put  in  the  system  in  a  horizontal  position.  The  pressure  drop  and 
the  energy  loss  were  observed .  The  drop  was  plotted  versus  the  flow 
rate  and  a  modified  loss  correlated  with  the  sediment  concentration.  An 
attempt  was  made  to  explain  the  results  theoretically.  It  is  necessary 
hereby  to  describe  the  resistance  of  the  grains  in  water  for  accelerated 


and  decelerated  motion . 
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II.  INTRODUCTION 

The  scope  of  this  report  was  to  develop  a  device  which  would 
measure  the  flow  rate  and  concentrations  in  a  solid-liquid  mixture . 

This  is  a  two-phase  problem,  with  water  as  the  liquid  phase  and  sand 
as  the  solid  phase . 

In  the  case  of  the  one -phase  flow  problem  a  sudden  change  in 
head,  caused  by  a  change  in  cross  section,  can  be  calibrated  to  indicate 
the  flowrate  in  the  system .  Application  of  this  phenomenon  is  made  in 
the  different  flowmeters ,  like  the  Venturi  meter ,  the  orifice  meter ,  etc . 

A  Venturi  meter  was  tested  for  its. applicability  in  a  two-phase 
system .  The  head  drop  should  give  us  some  information  about  the  flow- 
rate  and  it  was  hoped  that  another  measurement  on  the  same  instrument 
such  as  the  head  loss  across  the  Venturi  meter  might  give  us  some  in¬ 
formation  about  the  sediment  concentration.  This  Venturi  meter  appli¬ 
cation  is  discussed  in  Chapter  V  of  this  report . 

In  addition  to  the  Venturi  meter  another  system  was  believed  to 
give  the  desired  information  namely,  flow  rate  and  sediment  concentra¬ 
tion.  This  system  is  called  a  "Loop"  system.  The  Loop  system  con¬ 
sists  of  two  vertical  pipes  of  equal  length  and  of  opposite  flow  direction. 
Head  readings  are  obtained  across  each  of  the  vertical  pipes .  H  the 
pipe  length  is  the  same  in  both  sections,  a  summation  of  the  head  loss 
in  both  vertical  pipes  will  give  us  information  about  the  flow  rate .  The 
concentration  might  be  found  by  obtaining  the  difference  of  the  two  head 
loss  values .  This  should  be  equal  to  the  static  pressure  of  the  submerged 
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weight  of  the  solid  phase  over  the  length  of  both  vertical  pipe  sections . 
The  Loop  system  is  discussed  in  Chapter  IV . 


m.  EXPERIMENTAL  APPARATUS 


A  new  piece  of  experimental  equipment  completely  selfcontained 
and  independent  from  other  equipment  in  the  Hydraulics  Laboratory  was 
constructed.  This  was  necessary  for  following  reason:  (a)  all  pumps 
would  be  damaged  by  pumping  solid  material,  (b)  the  entire  water  supply 
would  become  contaminated  with  sediment,  and  (c)  the  system  would  need 
more  maintainance ,  since  it  would  require  frequent  cleaning. 

This  system  consisted  of  a  centrifugal  drive -pump,  a  driveline, 
a  jet  punip  for  the  sand -water  mixture,  a  dischargeline,  a  weighing  de¬ 
vice,  a  hydrocyclone  or  sand  separator ,. and  a  return  channel.  Each 
component  is  discussed  in  more  detail  below .  The  entire  system  is 
shown  in  Fig.,  III-l  and  Fig.nil.-2.  Fig.  III-l  is  a  general  photograph 
of  the  equipment  and  Fig.  III. -2  is  a  schematic  drawing  of  the  equipment. 
Drive  Pump.  The  drive-pump  was  located  one  floor  below  the  apparatus. 
This  was  necessary  to  facilitate  priming  this  pump .  Water  was  supplied 
from  the  return -channel  and  then  pumped  into  the  driveline .  The  pump 
was  driven  by  a  20  HP  440V  3  phase  motor  and  developed  a  capacity  of 
250  gpm  at  a  head  of  180  feet.  The  pump  remained  in  good  condition  at 
the  end  of  the  experiments,  indicating  an  efficient  separation  of  drive - 
water  and  sand  by  the  hydrocylohe , 

Drive  Line .  The  driveline  leading  from  the  discharge  side  of  the  drive - 
pump  with  a  nominal  diameter  of  2  inches .  A  gatevalve  was  installed  in 
this  line  and  was  used  to  regulate  the  discl^rge  from  the  jet  pump .  This 
valve  was  completely  open  for  maximum  velocity  in  the  discharge  line . 
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FIGURE  in -!"  EXPERIMENTAL  EQUIPMENT 

Driveline  Sediment-line  Loop  System 

Valve  Hydrocyclone  Manometers 

.  Venturi  ...  _  . 

<;u/i+rhee  fnr  Monometers  Weighing  Device 

centrifugal  Pump  Return  Channel 


m- 2  -  SCHEMATIC  OF  THE  EXPERIMENTAL  EQUIPMENT 
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The  jet  pump  was  provided  with  two  supply  lines  of  1-1/4  inch  each. 
Therefore,  the  2  inch  drive  line  was  divided  into  two  1-1/4  inch  drive- 
lines  near  the  jet  pump . 

Jetpump .  The  jetpump  was  developed  in  a  previous  study  for  pumping 
sand -water  mixtures .  The  main  flow  of  sand  and  water  passed  in  the 
pump  through  a  straight  2  inch  pipe  to  reduce  wear.  Three  jets  entered 
this  pipe  at  an  angle  of  approximately  lo®  from  the  outside  making  the 
main  pipe  the  mixing  chamber .  This  chamber  was  followed  by  a  dif- 
fusor  which  increased  the  main  pipe  to  a  3  inch  diameter.  By  increas¬ 
ing  or  decreasing  the  pressure  in  the  jets,  one  is  able  to  regulate  the 
main  flow . 

The  nozzle  ring  failed  once,  due  to  wear  near  the  nozzles .  A 
second  nozzle  ring  lasted  until  the  end  of  the  experiment .  The  jetpump 
is  described  in  Ref.  1 . 

Discharge  Line.  A  3  inch  discharge  line  leaving  the  diffusor  extended 
to  the  hydrocyclbne.  A  number  of  elbows  were  necessary  to  obtain  the 
desired  alignment  of  the  main  pipe  line ,  The  discharge  line  contained 
the  experimental  units  which  were  studied  under  two -phase  flow. 

Sampling  Device .  A  3 -way  valve  was  installed  in  the  discharge  line  and 
permitted  the  flow  to  be  interrupted  and  to  be  diverted  into  a  sampling 
device .  The  sampling  valve  was  designed  from  a  standard  pipe  -  T  and 
a  plunger ,  with  rubber  disks  sealing  against  pipe  surfaces  .  The  rubber 
disks  had  to  be  replaced  twice .  Attached  to  the  plunger  was  a  micro¬ 
switch  which  was  connected  to  a  timer  with  I/IO  of  a  second  least  count. 
The  sampling  discharge  was  diverted  into  a  weighing  tank  (94"  x  34"  x20"). 
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The  mixture  passed  through  a  box  with  a  screen  bottom  where  the  sand 
accumulated.  After  the  weight  of  the  mixture  was  obtained,  the  water 
was  drained  into  the  returnchannel  and  then  the  accumulated  wet  sand  was 
weighed . 

Hydrocyclon  or  Separation  Tank.  The  hydrocyclon  or  separation  tank 
making  use  of  centrifugal  force  separated  the  solid  phase  from  the  liquid 
phase .  Since  the  flow  in  such  a  device  is  rather  complicated  a  small 
model  was  first  constructed .  After  a  number  of  trials ,  a  satisfactory 
design  was  believed  to  be  found.  The  model  (Fig.  III-3)  was  operated 
by  a  positive  displacement  pump,  a  sand-water  mixture  supply  and  an 
emergency  jug,  to  catch  the  sand  passing  the  cyclone.  The  model  of  the 
hydrocyclone,  shown  in  Fig.  III-4,and  Fig.  III-5  was  made  of  lucite  with 
a  diameter  of  5  inches,  with  the  intake  pipe  at  2  inches  above  the  bottom 
and  the  outlet  about  5~l/2  inches  above  the  bottom .  A  metal  cone,  with  a 
diameter  of  3-1/2  inches  was  then  put  into  the  me  lei,  with  its  highest 
point  2-1/2  inches  above  the  bottom.  Beneath  the  hood  was  ample  space 
for  sediment  accumulation.  At  an  elevation  of  4-1/2  inches  a  baffle  ring 
was  installed .  This  ring  was  1/2  inch  wide .  The  wall  of  the  container 
was  kept  dirty,  for  increasing  the  wall  friction.  The  supply  bottle  was 
filled  with  sand  and  water .  By  means  of  adjusting  the  pipe  relative  to 
the  sediment  the  concentration  could  be  regulated.  The  tests  gave  the 
following  results .  The  critical  sand  size  with  no  accumulation  was  ob¬ 
served  to  be  d  =  0.32  mm .  The  test  was  started  with  a  sand  of  size 
d  =  1 .0  mm .  The  discharge  at  the  positive  displacement  pump  was  kept 
constant,  to  give  a  liquid  velocity  of  vl  =  1 .77  m/sec. 
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FIGURE  ni-3-  MODEL  EQUIPMENT 
Pump  and  Supply  Model  Emergency  Jug 

Motor  Bottle  Hydrocyclone 


FIGURE  m-4-  MODEL  OF  HYDROCYCLONE 
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FIGURE  in-5  -  MODEL  OF  SEPARATION  TANK 
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Assuming  that  Froude’s  law  governs  this  problem,  i.e.,  that 

viscosity  forces  are  of  minor  influence  only,  and  taking  the  velocity 

VL  =  1 .77  m/sec .  and  a  grain  diameter  of  d  =  3.2  x  10"4  the  Froude 
2 

number  vl  /gd  equals  980.  With  a  grain  diameter  of  d  =  1 .2  x  lO"^  we 
obtain  from  a  Froude  number  of  980  a  velocity  of  vl  =  3.43  m/sec. 
Therefore,  the  velocity  where  the  given  sand  grain  will  be  deposited  in 
the  hopper  of  the  separation  tank. should  be  smaller  than  3.43  m/sec. 

This  value  is  slightly  smaller  than  the  maximum  velocity  used  in  the 
prototype,  which  is  v^ax  =  3.74  m/sec.  The  operation  of  the  prototype 
proved  satisfactory .  The  intake  pipe  was  installed  so  that  the  flow  entered 
tangentially .  The  jet  circled  the  tank  and  then  the  water  rose  .  This 
rising  motion  continued  until  it  reached  the  tangential  outflow  opening  at 
the  top  of  the  separation  tank.  The  inner  and  outer  wall  of  the  hydro¬ 
cyclone  should  be  of  rough  texture  to  increase  friction  and  so  cause  the 
sand  to  settle .  The  sand  has  a  chance  to  accumulate  in  the  hopper  on 
the  bottom  of  the  hydrocycline .  The  motion  of  the  water  might  be.  ex¬ 
plained  as  a  secondary  current .  An  additional  baffle  ring  was  installed 
in  the  hydrocyclone  to  increase  the  separation  efficiency. 

The  prototype  was  constructed  as  shown  in  Fig.  III-6.  It  con¬ 
sisted  of  two  parts  such  that  the  upper  section  could  be  lifted  from  the 
lower  section.  To  the  upper  part  an  8  inch  pipe  was  welded,  which  con¬ 
tained  the  two  drivepipes  and  the  discharge  pipe.  The  diameter  of  the 
cyclon  was  2  feet  and  the  total  height  was  75  1/2  inches .  An  opening  on 
the  top  of  the  hydrocyclone  permitted  solids  to  be  added  and  the  valve 
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Outflow 

Pipe 


Inflow 

Pipe 


FIGURE  in-6- SEPARATION  TANK 
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on  the  bottom  of  the  hp|jper  permitted  the  system  to  be  drained .  The 
sediment  collects  in  this  hopper.  The  vertical  distance  between  the  in¬ 
take  and  the  outlet  pipe  is  26  inches .  The  inner  wall  of  the  hydrocyclone 
was  extremely  rough  as  formed  by  two  baffle  rings .  During  operation 
the  hydrocyclbne  was  open,  to  the  atmosphere. and  a  parabolic  water- 
surface  at  a  height  of  the  outlet  pipe  was  observed.  From  the  outlet 
pipe  the  water  flows  into  the  return  channel . 

Return  Channel .  The  return  channel ,  3  feet  wide  and  30  feet  long,  was 
subdivided  into  4  unequal  settling  basins .  The  first  one  was  the  longest 
and  the  last  the  shortest.  Thus  only  the  very  finest  particles  of  the  dust 
sizes  which  have  no  abrasion  effect  reached  the  centrifugal  pump .  The 
efficiency  of  this  settling  basin  therefore  was  satisfactory;  however ,  the 
whole  system ,  as  well  as  the  return  channel ,  has  to  be  cleaned  after  40 
hours  of  operation. 

Manometer  and  Pressure  Taps .  Standard  air -water  manometers  were 
used  to  measure  differential  heads .  The  pressure  taps  in  the  flow  lines 
were  made  as  follows;  A  1/16  inch  hole  was  drilled  through  the  pipe. 

A  standard  1/8  inch  coupling  then  was  welded  on  the  outside  of  the  pipe. 

A  1/16  inch  hole  was  sufficient  to  prevent  solids  from  entering  the 
manometer.  The  connecting  line  from  the  pressure  taps  to  the  manom¬ 
eter  was  saran  tubing . 

Water .  The  water  used  was  taken  from  the  general  water-supply  system 
in  the  Hydraulics  Laboratory .  At  the  beginning  the  water  was  clean  and 
clear.  After  operation  with  sand,  the  water  became  increasingly  turbid; 
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therefore,  it  was  necessary  to  clean  the  entire  system  frequently.  The 
water  temperature  was  recorded  and  found  to  vary  from  78®F  to  85*F 
during  the  entire  series  of  the  experiment. 


Si 
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IV.  THE  L(X)P  SYSTEM 

The  Loop.  System  consists  essentially  of  two  identical  vertical 
pipe  sections  with  opposite  flow  directions.  The  head  differences, were 
obtained  over  the  vertical  pipe  sections.  It  was  expected  that  the  sum¬ 
mation  of  the  head  differences  would  determine  the  flow  rate  while  their 
difference  would  determine  the  specific  weight  and,  therefore,  the  con¬ 
centration,  of  two -phase  flow  in  the  vertical  pipes.  The  experimeital 
verification  of  the  theory  is  presented  below . 

1 .  Theoretical  Considerations 

The  different  hydraulic  equations  for  the  liquid  phase  }ti  derived 
in  the  conventional  manner  as  follows,  and  the  effect  of  the  solid  phase  is 
then  considered  as  .an  external  force. 

1.1.  Basic  Equations.  .  . 

..  Fir;st  the."  d:q.ua'tion  .willlbeca:et:up-ih.a..yer.y  .general.  Wayi;  and  Ithen  we 
shall  cpnsider  them  in  an  application  for  the  Loop  system . 

1.11.  Continuity  Equation.  The  continuity  equation  states  that 
the  average  velocity  over  a  cross  section  multiplied  by  the  area  of  the 
cross  section  is  a  constant.  This  constant  is  equal  to  the  flow  rate  in 
volume  per  unit  time . 

(a)  First  this  equation  is  applied  to  the  solid  phase.  The  fraction 
of  area  which  is  occupied  by  the  solid  phase  is  written  as  A.c.  The 
velocity  of  the  solid  phase  may  be  desigmted  as  Vg  and  the  flow  rate  may 
be  denoted  as  Qs .  And  so  the  continuity  equation  reads : 


Qs  =  A.  o  •  Vs  =  const. 


(1) 


(b)  If  A*(l-c)  represents  the  area  of  the  liquid  phase,  vl  its 
velocity,  and  Ql  its  flowrate,  then  the  continuity  equation  for  the  liquid 
phase  is; 

Ql  =  A«(l-c)  •  VL  =  const.  (2) 

(c)  Also  the  total  flow  rate  of  both  phases  will  obey  the  law  of 
continuity.  That  is: 

Q^g  +  Ql  =  Q  =  const.  (3) 

Furthermore,  we  shall  establish  a  ratio  of  solid  flow  rate  to  total  flow 
rate  and  shall  call  this  ratio  at  . 

Qs/Q  =  a  (4) 

or  rewriting  eqviation  (4); 

(Qg  +  Ql)  01  =  Qflr  =  Qg  (4l) 

The  same  stated  for  the  ratio  of  liquid  flow  rate  to  total  flow  rate 


Ql/Q  =  (!“«) 

(5) 

and 

(Qs  +  Ql^i  “  =  Q(i  =  Ql 

(5I) 

Qg  and  Ql  in  equation  (4)  and  equation  (5)  can  be  replaced  by  equation  (1)  and 
(2) .  Such  combinations  result  in  the  following  relationships  . 


Qg  =  c .  A .  Vs  =  Q  .  O' 

(6) 

and 

Ql  =  (1  -c)  •  A  •  Vl  =  Q  •  (1-a) 

(7) 

We  may  sa^ 

{  that  equation  (3)  could  be  written  as 

Q  =  vA 

(3) 

Where  A  is  the  total  area  of  the  pipe  and  v  is  an  average  velocity  over 
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the  cross  section  A.  Combining  now  equation  (3  )  and  equation  (6) 
amounts  to: 
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V  = 


Q 

T 


(8) 


and  a  combination  of  equation  (3  )  with  equation  (7)  gives  us : 


A 


(1-c) 

(l-a) 


VL 


(9) 


For  all  practical  purposes  in  the  Loop  system  one  may  assume  that 
Vg  S  =  V  and  therefore  c  may  be  written  also  in  place  of  a.  Also, 
this  assumption  will  imply  uniform  concentration  c  in  the  entire  Loop 
system . 

1.12.  Momentum  Equation.  Newton’s  second  law  of  motion  states :  The 
time  rate  of  change  of  the  linear  momentum  of  a  body  is  directly  propor¬ 
tional  to  the  force  acting  on  the  body .  Suppose  m  is  the  mass  of  the 
particle,  v  the  instantaneous  linear  velocity  and  EF  is  the  sum  of  all  ex¬ 
ternal  forces  acting  on  the  particle,  this  equation  may  be  written  as: 

EF  =  -1  (m  .  v)  (10) 

dt 


Since  (m  •  v)  is  a  momentum,  this  equation  is  also  referred  to  as  the 
momentum  equation.  Assuming  there  is  no  change  of  mass  with  respect 
to  time ,  one  may  write : 


EF 


(11) 


The  problem  to  be  investigated  is  a  steady  flow  problem ,  so  that 

^  +  V  ~  ~ 

dt  ~dt  ds  ds 


% 


(12) 
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Using  the  result  of  equation  (12)  and  substituting  it  into  equation 
(11)  we  obtain; 

EF  =  m  .  (v  ~)  (13) 

ds 

The  application  of  this  equation  to  pipe  system  shall  now  be  investigated 
on  hand  of  Fig.  IV-1 . 


Area  A 


A  control  volume  element  is  considered  over  the  length  d;S.  The  cross 

sectional  area  will  be  denoted  i^ith  A.  The  liquid  shall  have  the  density 

Pl  g  will  be  the  acceleration  due  to  gravity.  The  mass  of  the 

liquid  involved  is  pL  •  A  .  ds  •  (1-c).  The  acceleration  according  to 

dvL 

equation  (12)  is  .  A  pressure  force  of  p  •  A  •  (1-c)  acts  at  the 

entrance  and  a  pressure  force  in  the  opposite  direction  of  (p+  ^  ds)  . 

A  •  (1  -  c)  acts  at  the  exit.  The  gravitational  force  of  pl  •  A  .  ds  . 
(1-c)  •  g  •  cos  0  acts  in  the  direction  of  ds .  The  frictional  forces 


♦* 


dFfr  will  oppose  the  flow .  The  particle  in  the  flow  may  not  have  the  same 
velocity  as  the  liquid.  They  will  then  exert  a  force  on  the  flow  which  may 
be  positive  if  it  retards  the  flow .  It  may  be  designated  as  dFp  .  We  shall 
now  enter  these  terms  into  equation  (11) .  So  that  we  have: 

Ps  .  ds  •  A  .  (1-c)  ]  .  VL  p.A.  -  j^p+  ^  ds]  .  A  •  (1-c)  + 

+  PL  •  g-  A.  ds  .  (1-c)  cos  0 -dF^j.  “dFp 

(14) 


The  equation  is  an  extended  dynamic  equation  with  friction.  Extended  by 
the  term  dFp.  This  is  an  external  force,  which  is  caused  by  the  suspended 
particles  in  the  flow .  Such  a  force  dFp  is  the  simplest  way  by  which  the 
effect  of  the  solid  phase  may  be  expressed . 

1 ,2  Application  of  the  Momentum  Equation  to  the  "LOOP"  System 

1 . 21  Upward  Flow  Section  or  Riser  Section .  Equation  (14)  shall 
be  applied  to  the  riser  section.  According  to  Fig.  IV-2  the  following  signs 
will  be  obtained  for  the  various  ternis  assuming  constant  concentration . 

0  =  (-dp  •  A~pl  ds*  A*  g)  ‘(l-c)  (15) 
-dFfr  -  dFp 

These  forces  are  the  net  pressure 
force,  the  gravity  force,  the  friction¬ 
al  force  at  the  pipe  wall  and  the  particle 
force.  The  change  in  velocity  dv,  due 
to  friction  and  particles  shall  be  small , 
so  that  the  left  side  of  equation  (14) 


Fig.  IV-2 


becomes  zero .  All  forces  in  equation 
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(15)  show  a  negative  sign,  that  is  to  say  they  are  opposing  the  flow. 
Integration  of  equation  (15)  gives 

O  =  (-Zip  •  A  -  PL  •  L  -  A*  g)  .(l-c)-Ffi.  -  Fp  (16) 

Here  L  is  the  length  of  riser  section .  C6iisidering^  Fp  which  represents 
the  weight  of  the  particles  under  water  ,this  gives : 

Fp  =  A- L-  c  .  (Ys  -  Yl)  (17) 

Using  equation  (IV)  and  substituting  it  in  equation  (16)  we  obtain 

O  =  (~z5p  •  A  -  PL  ‘A*  g)‘(l-c)-(p^^g.  A*L-c  -  Ffr;  (18) 

This  equation  was  derived  for  a  condition  of  upward  flow.  This  section 
of  the  Loop  system  is  called  the  riser  section.  The  other  part  of  the 
system,  the  downward  flow  section  or  downcomer  section,  is  next  dis- 
cvussed. 


1 .22  Downv/ard  Flow  Section  or  Downcomer  Section. 


m 
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1.23  Summary  of  the  Important  Equations  and  a  Logical  Design 
of  an  Experiment  .  The  most  important  equations  of  the  previous  chap¬ 
ters  are  equation  (18)  and  equation  (19)  which  for  the  riser  section  is 

Ffr  =  (-iipA-pL  •  L*A*  g)(l-c)  +  (Pg-pjJ  L*A*g‘c;  (18) 

for'  downcomer  section  is 

Ff  r  =  (-ApA  +  PL  L*A*g)(l-c)  +  (pg  -  pl)  L*A -g-c  •  (19) 


Dividing  both  equations  by  A*  y-^  and  we  obtain  for  the  riser  section: 

( 

TL 


A-Tl 


for  the  downcomer  section 


YL 


(20) 


!!!£_=  (-^  +L)(i  .c)  +  c  Tg-i r-L  L 
A-YL  .  Yl  Yl 


(21) 


The  riser  and  downcomer  sections  are  equal  in  geometry  and  flow  con¬ 
ditions  .  Corresponding  terms  in  equation  (20)  and  equation,  (21)  are  thus 
of  equal  value .  The  influence  of  the  various  terms  can  thus  be  separated 
by  mere  addition  or  substration  of  the  two  equations .  As  shall  be  shown 
numerically,  (this  chapter  part  3.3)  the  concentrations  in  equation  (20) 
and  equation  (21)  for  our  purpose, are. assumed  to  be  the  same.  The  first 
term  in  equation  (20)  and  equation  (21)  may  be  recognized  as  a  head 
across  the  length  L.  This  head  difference  will  be  obtained  across  the 
two  vertical  sections  at- the  Loop  system. 

1 .24  Head  Differences.  We  shall  first  write  down  the  equation 


for  the  head  difference  in  a  vertical  section.  Regardless  whether  the 
flow  goes  upward  or  downward  the  head  loss  is  due  to  the  pipe  roughness . 


We  may  write: 
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^  +  L  =  R  (22) 

YL 


where  R  stands  for  reading  on  a  manometer .  For  our  purpose  equation 
(22)  is  not  sufficient,  since  it  does  not  consider  the  suspended  solids  . 
We  can.  however,  modify  equation  (22)  such  that  we  consider  the  pres¬ 
sure  difference  due  to  the  solid  phase .  This  effect  of  the  suspended 
solids  act  with  a  positive  sign  in  the  riser  section  and  with  a  negative 
sign  in  tiie  downcomer  section.  Therefore,  the  modified  equation  (22) 
is 

(a)  for  the  riser  section: 


(-■^  -  L)(l-c)  -  L  (l£-l2k)  c  =  R„  (23) 

7L  YL  ^ 


(b)  for  downcomer  section: 


^  +  L)(l-c)  +  L  (  I^IIL  )  c  =  R^  (24) 

7l  ^ 

The  term  (  -  +  L)(l~c)  considers  the  frictional  effect  and  its 

T  “ 

magnitude  is  the  same  in  both  sections  and  is  termed  Rfj. .  The  other 

term  L  ^  )  c  is  due  the  suspended  solids,  and  is  termed  Hg^sp* 
7l 

So  equation  (28)  and  equation  (29)  can  be  written  as 
(a)  for  the  riser  section: 


^fr  ^susp .  ~  ^R 
(b)  for  the  downcomer  section: 

^fr  “  ^susp  ~ 


(25) 


(26) 


These  head  differences,  are  then  .substituted  into  equation  (20)  and  equa 
tion  (21)  and  we  get  for  the  riser: 
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hr  =R 
A.7l  ^ 


"  %r  ^susp 


for  the  downcomer: 


Ff*. 

7 —  ~  ~  %r  "  ®susp 

A-Tl 


Addition  of  the  two  equations  gives 


2 


Ffr 

A‘Yl 


2R. 


fr 


(27) 


(28) 


(29) 


The  effect  of  the  solid  phase  disappears  if  equation  (27)  and  equation 
(28)  are  multiplied  by  ILl. 

ym 


2  Jk.  =  (Rj^  +  Rd)  It  ( 30) 

A*  Tm  ‘7m 


where  ym  is  the  specific  weight  of  the  mixture  and  Rjj  7^ 

7m 


the. readings  in  head  of  mixture. 

Subtraction  of  equation  (27)  and  eqiuition  (28)  gives 


and  Rp  Tl 
7m 


0  =  R|^  ~  Rjj  =  3Rgjjsp  (31) 

Substituting  Rg^gp  for  L(  IS— It  )c  gives 

-  Rp  =  2  L  (-^~It  )  c  (32) 


With  equation  (30)  we  can  express  the  flow  rate  and  with  eqiiation  (32) 
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we  can  express  the  concentration. 

1.25' Summary.  A  mere  addition. of  the  manometer  readings  in 
head  of  mixture  gives  information  on  the  friction  loss  across  with  verti¬ 
cal  sections .  The  frictional  loss  in  Itself  provides  information  on  the 
flow  rate  of  the  mixture . 

The  subtraction  of  the  readings  of  the  manometers  indicates  a 
linear  increase  in  concentration.  Equation  (32)  states,  that  the  manoi- 
metric  difference  between  two  elevations  is  equal  to  the  immersed 
weight  of  the  material  held  in  suspension  between  these  two  elevations. 
Equation  (30)  and  Equation  (32)  are  compared  with'  experimental  data  in 
part  3  of  this  chapter . 

2.  Experimental  Datr. 

In  Chapter  III  the  experimental  equipment  and  procedure  was 
discussed.  The  sand  used  in  the  tests  are  described  as  follows: 

2 .1  Properties  of  the  Solid  Phase 

2.11.  3and  No .  2 .  In  most  experiments  a  Monterey  sand  Lapis 
Lustre  No.  2  was  used.  For  determining  the  specific  weight,  the  formu¬ 
la  jg  =  Wg/Vg  was  used.  For  seven  samples  an  overdry  weight  Wg  was 
determined .  The  volume  Vg  was  measured  in  a  calibrated  cylinder 
filled  partially  with  water ,  by  taking  the  displaced  liquid  as  Vg .  The 
ratio  of  Wg/Vs  was  fovtnd  to  be  2.607  kg/L.  The  ratio  of  wet  wei^t 
to^lry  weight  was  found  to  be  1,15.  A  sieve. analysis  (Fig.  IV -4)  showed 
that  dgQ  was  1 .15  mm  and  that  the  sand  was  fairly  uniform. 


I 

I  24 


M900tA 


FIGURE  IV- 4  -  MECHANICAL  ANALYSIS  OF 
SAND  N0.2  a  SAND  NO  H 
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2.12.  Sand  No ,  II .  This  sand  was  only  used  in  a  few  experiments  . 
The  specific  weight  was  determined  by  the  previous  procedure  and  was 
found  to  be  2.726  kg/L .  The  ratio  of  wet  weight  to  dry  weight  was  1.04. 
The  sieve  analyses  shown  in  Fig.  IV-4  gave  dgg  =  1.7  mm. 

2.13.  Summary  of  Sand  Characteristics  .  The  characteristics 

of  the  sands  used  in  the  tests  are  summarized  in  the  following  Table  IV- 1 . 


Sand 

Ys  kg/L 

dgo  mm 

wet  weight 
dry  weight 

Vo  4.4.1  «l/S 

Settl 

No.  2 

2.607 

1.15 

1.15 

0.11 

No.  II 

2.726 

1.70 

1 .04 

:  0.16 

Table  IVr  1 


2.2  Experimental  Procedure. 

The  experimental  equipment  was  previously  described  in, Chap¬ 
ter  III .  The  discharge  line  was  arranged  in  a  Loop  form .  This  means 
that  after  a  20  inch  horizontal  section  there  was  an  elbow  and  then  a  90 
inch  vertical  riser  section.  After  another  elbow,  a  short  52  inch  hori¬ 
zontal  section,  and  then  it  dropped  again  in  a  90  inch  vertical  section  to 
the  original  elevation.  -Figures  III^l  and  III-2  show  views  of  the  equip¬ 
ment. 

The  head  difference  across  the  vertical  sections  was  desired. 

It  was  obtained  in  both  cases  over  a  section  of  1 .5  m .  The  pressure 
taps  were  here  as  high  as  possible  in  the  riser  to  avoid  elbow  effects  as 
much  as  possible.  In  the  second  vertical  section  the  pressure  taps  were 
at  an  elevation  as  low  as  possible .  The  pressure  measured  at  any  lo- 
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cation  was  an  average  obtained  from  two  opposite  pressure  taps. 

2.21  Regulation  of  Flowrate .  The  discharge  in  the  test  section 
was  regulated  by  a  valve  in  the  jet  pump  drive  line  .  Thus,  one  could  ob¬ 
tain  various  pressures  in  the  jet  chambers.  This  was  the  most  effec¬ 
tive  device  to  change  the  flow  rate  in  the  experimental  pipe  without  wear 
in  the  valve.  As  one  might  expect,  also  by  increasing  or  decreasing  the 
concentration  the  flowrate  was  slightly  affected  at  a  given  setting  of  the 
valve . 

2.22  Regulation  of  Concentration.  The  concentration, in  the  ex¬ 
perimental  pipe  was  regulated  by  increasing  or  decreasing  the  total 
sediment  content  of  the  system . 

2.23,  Manometer  Readings .  The  head  across  both  vertical  pipe 
sections  was  obtained  by  two  40  inch  long  air-wate.’'  manometers .  Every 
set  of  runs  was  started  with  clear  water,  i.  e.:  zero  concentration. 

The  reading  on  both  manometers  were  the  same ,  as  expected .  Then 
gradually  the  concentration  was  increased  and  it  was  noticed,  that  at  the 
riser  section  the  reading  increased  and  at  the  downcomer  section  the 
reading  decreased  and  in  some  instances  gave  negative  values .  This 
behavior  is  in  accordance  with  the  theoretical  expectation  as  discussed 
in  part  1  of  this  chapter . 

The  manometer  scales  were  graduated  in  1/100  of  a  foot,  the 
1/1000  ft  was  estimated.  Minor  fluctuations  always  existed.  They 
were  damped  by  friction  in  the  saran  tubing  and  by  keeping  the  petcocks 
on  the  manometers  only  partially  open.  At  some  intermediate  concen- 
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trations ,  and  at  certain  flow  rates  the  damping  was  not  sufficient  and 
an  average  had  to  be  estimated .  It  is  believed  that  this  is  a  problem  of 
instability  when  the  solid  phase  probably  moved  in  clouds . 

The  manometers  were  read  before  and  after  samples  of  the  mix¬ 
ture  were  taken.  The  sampling  procedure  is  explained  in  section  2.24. 

The  difference  in  the  manometer  readings  before  and  after  sampling  was 
small,  since  the  sampling  time  was  short.  The  averages  of  the  read¬ 
ings  were  then  determined  for  each  vertical  section. 

The  readings  for  the  riser  and  downcomer  section  are  tabulated 
in  column  4  and  5  of  Appendix  I . 

2 .24  Sampling  Procedure .  The  sampling  of  sediment  was  ac¬ 
complished  after  equilibrium  was  reached,  i.e.,  when  the  concentration 
remained  constant.  This  could  be  checked  on  the  manometer  readings.. 

The  sample  was  taken  from  the  flow  by  means  of  a  three-way  valve, 
which  diverted  the  flow  for  a  desired  time  into  a  weighing  tank ,  The 
weighing  was  done  as  described  in  Chapter  III.  The  sampling  time  used 
in  the  experiment  ranged  from  0.8  seconds  to  2.2  seconds.  The  reason 
for  the  short  sampling  time,  was  to  prevent  an  excessive  decrease  of 
concentration  in  the  system.  The  weight  in  kg  of  mixture,  of  sand  and 
of  water  then  was  recorded.  The  weight  in  kg  of  mixture  ranged  from 
10  kg  to  17  kg -and  the  weight  in  kg  of  sand  from  0  kg  to  6  .5  kg. 

2.3  Determination  of  Flow  Rate  and  Concentration 

2.31  Determination  of  Flow  Rate;  The  flow  rates  were  deter  - 
mined  from  the  sampling  weights  and  times  .  Since  the  weight  of  sand 
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was  obtained  as  a  wet  weight,  the  values  had  to  be  corrected  to  dry 
weight.  The  correction  factor  is  given  in  Table  IV.  1.  In  the  same 
table  Yrn  is  given  for  both  sands .  Dividing  the  dry  weight  by  y  results 
in  the  volume .  The  flowrate  was  obtained  by  dividing  the  volume  by  the 
time,  and  then  presented  in  liter /sec.  The  flow  rates. are  given  in 
column  1  of  Appendix  1 . 

2.32  Determination  of  Concentration;  We  want  to  have  the 
concentration  in  volume  percentage.  This  was  obtained  as  the  ratio  be¬ 
tween  sediment  flow  rate  and  total  mixture  flow  rate  .  This  was  done 
for  each  run  and  is  tabulated  in  column  7  of  Appendix  1 . 

2.4  Summary  of  Experiment 

From  our  experiment  we  obtained  the  readings  at  an  air -water 
manometer  in  head  of  water  for  the  riser  and  the  downcomer  section. 
Furthermore ,  the  flowrate  of  the  mixture  was  obtained  in  liters  per 
seconds  and  the  concentration  in  volume  percent .  These  four  values 
are  tabulated  for  each  run  in  Appendix  1 . 

3 .  Comparison  of  Experiment  and  Theory 

3 .1  The  Concentration 

We  shall  now  compare  the  experimental  results  with  equation 
(32).  Rewriting  this  equation  gives 

Rr  -  RD  =  2'L.  )  .  c  (32) 

The  value  (2  L)  is  3  meters  .  Rr  and  Rjj  are  obtained  as  manometer 
readings  in  head  of  liquid . 

For  the  two  different  sands  (  — — ^ 

rL 


)  was  1 .607  for  sand  No .  2 
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and  1  .726  for  sand  No.  II.  Using  these  experimental  values  we  obtain 
for  sand  No .  2 

%  ~  Rj)  =  4.821  o  (33) 

for  sand  No.  n 

Rjj  -  Rd  =  5.178  c  (34) 

Both  eqimtions  plot  as  straight  lines  for  -  Rj-^  versus  conpentration. 
Equation  (26)  is  plotted  in  Fig.  IV-5  .  Also  in  this  figure  are  plotted 
about  350  experimental  values  for  sand  No,  2 .  It  is  apparent  that  the 
experimental  data  fit  the  theoretically  line  quite  well  despite  of  the  rather 
large  scatter.  This  is  especially  true  at  the  lower  and  intermediate  con¬ 
centrations,  but  at  higher  concentrations  (>  12%)  the  experimental  data  niay 
deviate. systematically  somewl^at  from  the  predicted  line.  No  explana¬ 
tion  is  given, for  this.  In  Fig.  IV-6,  the  experimental  runs  with  sand 
No.  U  are  compared  with  equation  (34).  Also  here  it  is  noticed,  that 
the  data  are  in  good  agreement  with  the  theoretically  predicted  line .  The 
experimental  data  lised  in  Fig.  IV-5  and  Fig.  IV-6  are  tabulated  in 
Appendix  1 ,  For  these  two  figures  it  was  necessary  to  calculate  the 
differences  of  the  two  manometer  readings,  at  the  riser  and  at  the  down¬ 
comer  section,  as  tabulated  in  column  6  of  the  Appendix  1 ,  starting  with 
zero  at  zero  concentration  and  increasing  with  concentration . 

3 . 2  The  Flow  Rate  .  The  theoretical  equation  derived  in  part  1 
of  this  chapter  is ; 

2  .  pi-  =  (Rr  +  Rn)  X  (30) 

Arm  "  “  tm 

The  right  side  of  this  equation  is  the  manometer  reading  in  head  of  mix- 


1,1^ 


A 

■ 

A''  ° 

A/*" 

0  * 

o  o 

<1 

D 

o 

0  □*,  ° 

^0 

*  % 

/  j  ® 

o  ♦  + 

*  0 
t3  . 

Sond  No.  2 

« 

d„*  1.15  mm 

+  t* 

0 

o  Experiment  60 

.  »  50 

.**o 

a  •  40 

o 

•  "  30 

0  ** 

)( 

♦  .  20 

E 

auction  of  heov 

V  line  Rr-Rb«  4,821c 

(uniform  concentrotlon  in  entire  cyttem) 

(chonge  in  conctntrotion  in  system) 

''imo!"  Rr-Rd  “♦•67  C 

<0f  Rr-Ro  •4.58  c 


Concentrotlon  c  -(Volume  %) 


0 

2  4  6 

8  10  12 

14  16  18  20  i 

HMOSC 

FIGURE  IV  -5 

-CONCENTRATION  vs 

MANOMETER  READING  DIFFERENCE 

LOOP  SYSTEM 

[Volume  x]  vt 

[Meters  or  fiMt  of  Water] 

iffercnc*  in  Manometer  Reading  in  meters  of  Water 


2.4 


.0.73 


ture .  This  value  of  (Rjj  +  Rj))  was  determined  for  each  run  and  is 
shown  in  column  3  of  Appendix  1  in  head  of  water  .  The  values  of  this 
term  were  mxdtiplied  by  ylym  and  tabulated  in  column  2,  Appendix  1 . 
The  specific  weight  of  the  mixture  was  obtained  as  follows ; 

Ym  =  YL  +  (rs  -  Yl)  c  (31) 

The  left  aide  of  equation  (30)  is  the  friction  loss  over  both  the  riser  and 

the  downcomer  sections .  The  frictional  loss  is -generally  speaking- 

proportional  to  the  square  of  the  velocity  and  flow  rate.  Therefore,  the 

right  side  of  equation  (30)  was  plotted  versus  the  flow  rate .  The  circled 

points  shown  in  Fig.  IV-7  were  obtained  for  clear  water.  All  other 

points  represent  various  concentrations .  They  show  some  scatter 

/ 

around  the  line  for  clear  water.  On  Fig.  IV-7  experimental  values  of 
both  sands  are  plotted.  It  is  expected  that  the  points  follow  on  a  log- 
log  paper  a  line  with  an  apprt  cimate  slope  of  2.  This  is  shown  on  Fig. 

IV-8,  where  the  slope  of  the  best  fitting  straight  line  of  Fig.  IV-7  was 

\ 

2.17. 

3 . 3  Confirmation  of  the  Assumption  of  Constant  Concentration 
in  Entire  Loop  System ; 

That  the  concentration  in  the  riser  and  in  the  downcomer  section 
are  equal  was  assumed  in  part  1.23  of  this  chapter  and  is  proved  as 
follows ;  by  applying  the  continuity  equation  for  the  solid  phase  in  the 
riser  and  the  downcomer  section.  Therefore,  the  same  amount  of  dis¬ 
charge  which  first  goes  up  in  the  riser  section  must  then  come  down  the 
downcomer  section .  We  write  therefore , 


FLOW  RATE  vs  MANOMETER  READING  SUMMATION 
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Qg,  in  riser  =  Qg,  in  downcomer  (35) 

For  Qg  we  use  equation:  (1)  which  combined  with  equation  (35)  gives 

(A*  c  •  Vs)in  ~  ^A*  c  •  vs)in  downcomer 

The  area  A  appears  on  both  sides  therefore  equation  (36)  reduqes  to 

(c  .  Vs)jj^  riser  ~  '  ^s^in  downcomer 

The  solid  phase  velocity  in  the  riser  section  can  be  expressed  according 
to  Fig .  IV- 2  as  being 

vs  in  riser  =  vl  -  Vgettl 

The  means  that  the  solid  phase  will  move  slower  than  the  liquid  phase 
by  an  amount  of  its  settling-velocity.  From  Fig,  IV-3  we  may  obtain 
that  in  the  fall  section  the  sand  velocity  will  be  larger  by  the  amount  of 
the  settling  velocity .  Then  we  write : 

Vg  in  downcomer^  =  vl  +  v  getli 

Combining  equation  (38)  and  equation  (39)  with  equation  (37)  we  obtain 

''^L  ■  %et'i:l 

- 7 -  =  (40) 

’'L  +  ’"solil  'R 

Where  c'j)  is  the  concentration  in  the  downcomer  section  and  the  concen¬ 
tration  in  the  riser  section  is  cjj.  The  influence  of  the  settling  velocity 
will  cause  the  deviation  of  cr  from  cq.  By  means  of  numerical  values, 
taken  from  the  experiments,  it  is  shown,  that  the  concentration  does  not 
change  significantly  in  the  entire  Loop  system.  The  following  table 
(Tab .  IV- 2)  illustrates  this  fact  for  the  two  sands  used  in  the  tests . 
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Table  IV-2 


Sand 
No.  2 

Vsettl  =  O.llm/s 

''Lmax  =  3.7m/s 

''Lmin 

c-'d/c.r  =  0.945 

,c.jj/cjj  =  0.909 

Sand 
No.  n 

''^'settl  =  0.16  m/s 

VLmax  =  3.7  m/s 

VLmin  =  2.2m/s 

cd/cR  =  0.920 

cj)/cR  =  0.866 

The  maximum  and  minimum  test  velocities  were  selected.  All  other 
values  will  be  found  between  this  two  extreme  velocity  values .  We  now 
replace  c  in  equation  (32)  by  (cr  +  cd)  .1/2  and  v/ill  investigate  how  far 
this  affects  equations  (26)^^hd  (27).  Also  this  is  summarized  in  Table 
IV-,3 . 


Table  IV -3 


Sand 

Velocity 

If  CD  =  CR 

If  Assumed  no  Change  in 
Concentration  in  the  Whole 
Loop,  i.e.  Cd  =  CR 

No.  2 

Vmax 

4.67  c 

4.821  c 

^min 

4.58  c 

(Eq.  33  ) 

No.  II 

^max 

4.95  c 

5.178  c 

Vmin 

4.83  c 

_ 

(Eq.  34) 

From  Table  IV-3  we  conclude,  that  at  the  small  liquid  velocities,  the 
deviation  from  the  eqioation  (33)  or  equation  (34)  are  the  largest,  yet  this 
deviation  is  quite  small.  It  can  be  observed  by  examination  of  Fig.  IV-4 
and  Fig.  IV-5,  that  our  assumption  of  a  uniform  concentration  all  through 
the  Loop  system  was  justified,  and  the  effect  of  the  settling  velocity  is 


of  insignificant  influence  in  this  investigation. 
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V.  .  VENTURI  METER 

In  chapter  IV  the  Loop  system  was  considered .  A  successful 
result  was  expected  and  the  expectation  was  satisfied  since  it  permitted 
the  determliiation  of  the  flowrate  and  the  concentration  by  manometer 
readings  only.  But  the  Loop  system,  due  to  its  rather  awkward  shape, 
might  not  be  the  best  for  practical  purposes.  .  For  this  reason  another 
design  was  investigated.  It  was  a  modified  Venturi  meter. 

1 ,  Experimental  Results 

The  applicability  of  a  Venturi  meter  in  a  two -phase  flow  problem 
is  considered.  .We. know  that  in. a  single-phase  flow  the  pressure  drop 
in  the  contracted  section  can  be  correlated  with  the  flow  rate.  ,  In  order 
to  obtain  some  additiomil  information  on  the  concentration  additional 
measurements  were  needed.  As  such  it  was  decided  to  obtain  in  the 
two-phase  flow  problem  the  head-loss  across  the  Venturi  meter;  how¬ 
ever  ,  it  is  necessary  to  explain  the  experiment  in  order  to  establish  a 

I 

theory . 

The  Venturi  meter  as  used  (see  Fig.  Ill- 1)  was  placed  in  a  hori¬ 
zontal  section  of  the  3  inch  pipe  line  system .  The  throat  diameter  of 
the  Venturi  was  2-1/8  inches,  A  longitudinal  section  through  the  Ven¬ 
turi  is  shown  in  Fig.  V-1 .  The  Venturi  meter  was  especially  paachined 
for  this  purpose  and  special  care  was  taken  to  obtain  a  smooth  transition 
from  the  pipe  to  the  Venturi .  The  pressure  drop  was  recorded  between 
tap  locations  (1)  and  (2)  by  means  of  an  air-water  manometer.  The 
head  difference,  converted  into  head  of  mixture,  was  then  plotted  hi 


SPACING  OF  THE  PRESSURE  TAPS 


% 


FIGURE  V-l  -  CROSS  SECTION  THROUGH  VENTURI  METER 
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graph  Fig.  V-2  and  2A  versus  the  flow  rate  and  is  also  tabulated  in 
Appendix  2 .  The  points  for  the  various  concentrations  fell  on  the  clear 
waterline  for  sand  No.  2,. as  may  be  seen.  For  Sand  No.  II  they  do  not 
seem  to  correlate  quite  as  well  with  the  clear  water  line,  especially  for 
the  high  concentrations .  In.  Fig.  V- 3  the  clear  water  line  is  replotted 
on  log-log  paper  and  seems  to  follow  quite  well  a  straight  line  with  the 
slope  2.26.  This  is  in  accordance  with  the  theory,  that  for  turbulent 
flow  any  head  difference  is  proportional  to  the  square  of  the  velocity,  as 
long  as  the  sediment  motion  deviates  only  insignificantly  from  that  of 
the  water . 

Another  reading  was  obtained,  which  indicated  the  energy -loss 
across  the  Venturi  meter.  The  pressure-taps  (a)  and  (b)  (see  Fig.  V-1) 
were  used  for  this  purpose .  To  avoid  any  local  effect  of  the  Venturi 
meter  the  pressure  taps  were  located  about  4"  away  from  the  beginning 
and  the  end  of  the  taper.  This  seems  to  be  quite  important  in  the  down¬ 
stream  section,  where  we  must  take  tlie  possibility  in  account  that  some 
separation  could  occur.  This  was  checked,  by  having  a  series  of  pres¬ 
sure  taps  along  the  exit  of  the  Venturi  meter .  After  study  of  the  results , 
obtained  with  different  concentrations,  it  was  decided  to  use  tap  (b) . 

The  head -loss  .across  the  meter  between  (a)  and  (b)  was  measured  by 
another  air -water  manometer .  The  reading  may  be  called  "b" .  It  was 
observed  that  "b"  increased  with  concentration.  For  this  reason  "b" 
was  plotted  against  concentration  and  the  flow  rates  noted  at  all  points 
(Fig.  V-4  and  4A) .  The  values  of  "b"  and  concentration  c  are  tabulated 
in  Appendix  2.  One  is  able  to  draw  through  the  different  experimental 
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points  lines  of  equal  flow  rate  or  equal  velocity.  All  of  these  equal-flow  lines, 
ranging’fr-qm  9  1/2  (.JseicAo  Ih  ifaec,  Have  the^ same-slope .  We:  now  assume 
that  the  losses  through  a  venturi  due  to  contraction,  expansion  and  wall 
friction  are  just  a  function  of  the  velocity.  This  part  of  the  loss  shall  be 
called  "bo"  .  "bo"  may  thus  be  determined  as  the  energy  loss  at  zero 
concentration  and  shall  be  subtracted  from  the  total  loss."b".  The  re¬ 
mainder  "b-bo"i  tabulated  in  Appendix  2,  was  plotted  in  Fig.  V-5  against 
concentration.  This  graph  shows  that  regardless  of  the  flow  rate  the 
value  "b-bo"  is  uniquely  related  to  the  concentration. 

We  shall  go  back  to  Fig.  V-4,and  4A  and  reason,  why  some  of 
the  points  do  not  follow  the  equal  flow  rate  lines  .  It  was  observed  at 
pressure  readings  just  upstream  of  V:enturi  meter  that  the  velocities 
increased  by  adding  more  solid  material.  This  means,  that  at  a  low 
velocity  the  water  is  unable  to  carry  all  of  the  solid  phase.  Part  of  it 
settles  and  cauF^s  a  decrease  in  cross  section,  which  may  be  recorded 
on  a  manometer  as  an  increase  in  velocity.  On  the  other  hand  the  ma¬ 
terial  which  accumulates  on  the  bottom  of  the  pipe  will  not  pass  through 
the  Venturi  meter  and  therefore  the  head  "b-bo"  is  expected  to  decrease, 
and  also  the  concentration  should  decrease.  The  first  effect,  a  decrease 
ih."b-bo"can  be  seen  clearljr  in.. Fig.  V-4  and  4A.  The  second  effect  a 
decrease  in  the  concentration  moved  through ,  does  not  show  in  the  graph 
of  Fig.  V-4  and  4A.  This  is  due  to  the  fact,  that  the  concentration 
plotted  is  the  fraction  of  solid  phase  involved  in  the  investigation.  By 
comparing  Fig.  V-4  and  Fig.  4A  we  notice,  that  in  Fig.  V-4A  accu- 
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mulation  occurs  at  lower  flow  rates .  The  reason  of  this  probably  is  due 
to  the  larger  grain  size . 

From  Fig.  V-5  we  can  obtain  the  relationship  between  "b-bo” 
and  the  concentration.  This  is 

(b-bo)  =  1.95  .  1^.  c  (41) 

Here  ly  is  the  length  over  which  "b-bo"  recorded  and  is  a  geometric 
constant  of  the  particular  instrument . 

c  is  the  concentration  and  1.95.  is  a  dimensionless  parameter. 

This  parameter  may  be  expected  to  change  with  the  type  of  sediment,  but 
the  experimental  results  indicated  the  same  parameter  for  both  kinds  of 
sand,  with  a  slight  variation  in  diameter  (dso  =1.15  mm  and  1 .7  mm)  and 
in  density  {y  =  2.60  kg/L  and  2.72  kg/L),  it  is  believed  that  this  parameter 
is  a  function  of  grain  geometry  as  well  as  of  the  density.  Further  re¬ 
search  must  develop  this  relationship .  In  addition  to  the  fact  that  the 
parameter  could  change  with  the  characteristics  of  the  material,  it  will 
also  change  with  geometry  of  the  Venturi  meter,  i.  e. ,  with  the  ratio 
Al/A2and£y, 

The  experiment  described  above  was  done  with  water  as  the  liquid 
phase  and  sand  No.  2  and  No.  II  as  the  solid  phase .  As  stated  in  Chap¬ 
ter  IV,  the  sand  No.  2  had  a  mean  diameter  of  1 .15  mm,  the  other  one 
had  a  mean  diameter  of  1 .7  mm.  It  should  also  be  mentioned,  that  the 
flow  rate  and  the  concentration  were  obtained  by  the  calibrated  Loop 
system,  which  was  arranged  in  series  with  the  Venturi.  This  arrange¬ 
ment  accelerated  the  process  of  obtaining  the  experimental  data. 


m 


FIGURE  V-5  -  VENTURI  METER 
b- bo^fTjeters)  vs  Concentration-(Vol.%) 
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2.  Theory 


2.1  The  i^essur'e  Drop 

We  shall  try  to  develop  a  theoretical  explanation,  why  the  pres- 
rues  drop  reading  caused  by  the  flow  of  a  mixture  in  a  Venturi  permits 
determination  of  the  flow  rate.  For  this  purpose  let  us  consider  Fig. 
V-Sci .  We  will  assume  that  the  Venturi  has  a  horizontal  axis  . and  that  the 


Fig.  V-6c< 


mixture  is  incompressible.  Fimthermore,  we  shall  assume  that  the 

concentration  will  not  vary  from  cross  section  Aj  to  cross  section  A2. 

We  shall  neglect  frictional  effects  and  assume  a  uniform  velocity  distri- 

« 

bution  at  the  inlet  and  at  the  throat.  In  such.a  case  the  energy  equation 
for  steady  flow  gives  the  relation  r 


H  +  iBl-  IL+  ^"^2 

Tm  2g  Km.  Tg 


(42) 


where  pi  and  p2  are  the  heads  in  points  (1)  and  (2),  Yjjj  is  the  specific 
weight  of  the  mixture  and  vm  is. taken  to  be  the  velocity  of  the  mixture-. 
In  addition  to  the  energy-equation  we  can  use  the  continuity  equa¬ 


tion,  which  reads; 
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Qm  ~  -^1  ~  Vni2  -^2 


(43) 


here  Qjji  is  the  total  flow  rate. 

Combining  equation  (42)  and  (43)  we  obtain 
Qm“  Cy 


A2 


2g(p2  -  Pi) 


(44) 


Ym 


Since  the  derivation  has  been  done  for  an  ideal  fluid,  the  factor  Cy  stands 
for  additional  real  fluid  effects  and  approaches  0.98  at  high  Reynolds 
numbers  .  Equation  (4)  could  be  rearranged  so  that  it  reads : 


Qm 


A2 


V  l-(Al/A2)2 


.  Y2g 


P2  "  PI 


Ym 


(45) 


The  bracket  term  on  the  right  side  of  this  equation  is  a  constant  for  our 

P2'  "  PI 

experiment ,  and  the  term  - ^  is  equal  to  the  head  of  mixture , 

Ym 

This  is  good  for  the  turbulent  flow  and  we  see  immediately  from  equation 
2 

(45)  that  Qm  ^  amjWhere  a^  is  the  pressure  drop  in  meters  of  mixture. 
As  a  summary  we  can  say  that  the  laws  which  explain  the  Ven¬ 
turi  effect  in  a  liquid  *  also  explain  the  flow  of  a  mixture  through  a  Ven¬ 
turi  meter, if  certain  assumptions  are  made.  The  pressure -drop  read¬ 
ing  has  to  be  taken  in  terms  of  the  mixture . 

The  experimental  results,  tabulated  in  Appendix  2.  and  plotted  in 
Fig.  V-2  and  2A  seem  to  prove  this.  The  slight  systematic  deviations 
in  Fig.  V-2A  indicate,  however,  that  the  concentration  as  well  as  the 
solid  particle  size  has  to  be  studied  more  extensively . 


.2 ...2  ..Energy  Loss 
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2.2.1  The  Drag  Equation  .  An  attempt  shall  be  made  to  explain 
the  energy  loss  across  the  Venturi  meter .  In  Fig.  V“4  and  4A  the  total 
energy  loss  is  plotted  against  concentration .  This  total  energy  loss  may 
be  called  "b".  It  includes  the  firctional  loss,  the  losses  due  to  contrac¬ 
tion  and  expansion  and  a  certain  loss  due  to  the  solids  .  Furthermore,  it 
was  shown. in  part  1  of  this  chapter,  that  under  the  assumption  that  the 
friction -expansion-  and  contraction, losses  are  linearly  superimposed 
on  the  concentration  effect,  these  two  effects  can  be  separated.  In  Fig. 
V-5  this  was  done  and  it  is  seen,  that  the  effect  due  to  the  concentration 
is  in  this  investigation  independent  of  velocity  or  flow  rate  with  the 
veloeitieB.’ranging  from.a  .;maximum:of  4:m/.sec  tO;  a  minimum 
of  2  m/sec .  The  relationship  between  the  head  and  the  concentration  c 
was  established  in  equation.  (41); 

(b-bo)  =  1.95.  Iv  •  c  (41) 

The  solid  particles  had  been  assumed  to  move  with  the  same  velocity  as 
the  liquid  phase  in  the  case  of  steady  flow.  This  is  to  say,  that  there 
exists  no  relative  motion  between  the  two  phases .  In  the  Venturi  meter 
we  have  significant  accelerations  with  the  corresponding  pressure  gradi'-.. 
ents.  Since  the  mass  of  a  solid  particle  is  bigger,  than  that  of  the  dis¬ 
placed  fluid,  its  .acceleration  in  a  pressure  gradient  will  be  smaller  than 
that  of  the  liquid.  This  will  cause  a  relative  motion  between  the  solid  and 
liquid  phase .  This  relative  motion  will  be  different  in  each  cross  section 
of  the  Venturi  meter.  A  certain  amount  of  energy  is  necessary  to  main- 


%lt 
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tain  this  relative  motion.  And  this  seems  to  be  the  additional  energy- 
loss  due  to  the  presence  of  the  solid  phase.  We  shall  express  the  re¬ 
sistance  of  Ihe  sand  in  water  with  a  drag  equation.  Then  we  get: 

(b-bo)  =  -SESBL  (46) 

r  •  A 

A  detailed  examination  of  this  drag,  equation. shows  that  the  relative  mo¬ 
tion  changes  throughout  the  Venturi  or  in. other  words  the  solid  particle 
is  accelerated  with  respect  to  the  liquid.  Obviously  the  steady-state 
drag  equation  is  insufficient  for  this  type  of  problem .  The  concept  of 
virtual  mass  was  established  in  the  18th  Century.  Later  Bessel  (Ref. 

11)  defined  the  virtual  mass  as  being  the  mass  of  fluid  displaced  by  a 
body  multiplied  by  a  factor  k.,  the  virtual  mass  coefficient.  In  Ref.  10 
Prandtl  and  Tietjens  explain  this  phenomenon  as  follows : 

"For  an  accelerated  motion,  it  was  seen  that  potential  flow  does 
lead  to  a  resistance .  In  order  to  accelerate  a  sphere  in  an  ideal  fluid 
it  is  not  only  necessary  to  exert  a  force  equal  to  the  product  of  the  mass 
of  the  sphere  and  its  acceleration,  but  an  additional  force  is  required  to 
accelerate  the  mass  of  the  fluid  particles  set  in  motion  by  it.  From  the 
above  equation  for  the  resistance,  it  is  seen  that  this  additional  force 
is  equal  to  the  product  of  the  acceleration  of  the  sphere  and  the  mass  of 
an  amount  of  fluid  of  half  its  voliune .  The  .apparent  increase  in  mass  of 
various  bodies  depends  on  the  shape  and  on  the  direction  of  motion .  For 
instance,  for  the  two-dimensional  flow  around  a  circular  cylinder,  the 
apparent  increase  in  mass  is  equal  to  the  full  mass  of  the  cylinder  in 
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liquid."  For  a  number  of  other  bodies  the  added  mass  coefficient  is 
summarized  in  the  paper  by  Iversen  &  Balent  (Ref.  2) .  These  added  mass 
coefficients  are  of  limited  importance  since  the  real  fluid  behaves  dif¬ 
ferently  than  the  ideal  fluid.  In  1924  Lunnon.(Ref.  3  and  Ref.  4)  reported 
from  his  own  experiment  that  this  coefficient  varied .  His  investigation 
with  spheres  showed  a  variation  from  k  =  0 .5  to  k  =  2 .0,  where  the  one 
obtained  with  potential  flow  was  reported  to  be  k  =  0.5 .  Iversen  & 

Balent  (Ref.  2)  reported  in  their  paper,  that  their  research  on  a  disk  re¬ 
sulted  in  an  added  mass  coefficient  of  ,k  =  0.637  to  k  =  5.0,  where  the 
lower  value  is  the  one  obtained  with  potential  theory.  The  same  authors 

were  able  to  obtain  by  dimensional  analyses. an  ^'acceleration  parameter" 
A  •  D 


and  established  fairly  successfully  a  relationship  between  the 
added  mass  coefficient  and  this  acceleration  parameter .  This  result  is 
shown  in  Fig.  V-9.  Also  shown  in  the  same  graph  are  Bugliarello 's 
(Ref.  5)  research  with  spheres.  The  variation  of  the  added  mass  co¬ 
efficient  was  from  k  =  0.4  to  k  =  100.  He  failed  to  find  any  relationship 
between  his  k  values  and  the  acceleration  parameter  .  With  the  few  re¬ 
sults  available  one  may  see  the  complexity  of  a  problem  dealing  with 
non-steady  flow . 

We  shall  try  to  set  up  an  equation  for  total  drag.  The  steady 
state  drag  equation  can  be  written  as  : 

DRAGi  =  Ki  Up2  (47) 


Additional  to  this  drag,  we  will  have  to  introduce  the  acceleration  drag. 
This  can  be  done  in  two  ways . 


u 
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The  first  way  might  be  as  follows:  To  the  steady  state  drag,  the 
acceleration  drag  has  to  be  added  and  we  get  : 

Drag2  =  Kj  Up^  +  k(VsPL)  Up  (48) 

For  a  given  body  geometry,  Ki,  in  equation  (47)  and  equation 
(48)  is  a  function  of  Reynolds  number  and  k  shall  be  a  function  of  the 
acceleration  parameter.  Vg  is  the  displaced  volume  of  liquid. 

'  The  second  way  to  take  in  account  the  acceleration  forces  is  to 
modify  the  steady  state  drag  equation.  In  this  case  equation  (47)  is: 

Drags -KsUp^  <^9) 

Here  K3  is  not  only  a  function  of  Reynolds  number,  but  also  of  the  ac¬ 
celeration  parameter . 

It  seems  reasonable  to  use  equation  (49)  where  the  acceleration 
effect  is. of  minor  importance. 

2.2.2  Equation  of  Motion .  It  was  seen  at  the  beginning  of  part 
2.2.1  of  this  chapter,  that  there  will  be  a  relative  motion  between  the 
solid  and  liquid  phase  throughout  the  Venturi  meter .  This  relative  ve¬ 
locity  is  responsible  for  the  existing  drag,  or  in  other  words  it  is  re¬ 
sponsible  for  the  additional  energy  loss,  which  is  termed  "b-bo".  This 
relative  velocity,  averaged  over  the  length  of  the  Venturi  meter  where 
the  loss  bbtaihe'diis  to  be  determined .  For  this  reason  the  differential 
equation  is  established  for  a  solid  particle  being  moved  by  the  pressure 
gradients  of  the  liquid  and  exerting  a  certain  drag.  The  differential  equa¬ 
tion  of  motion  is  in  such  a  case 

•  • 

pg.  Vs-  Vg  =  PL-  Vs-vL  -  Drag 


(50) 
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For  the  drag  equation  (49)  is  introduced  and  v  is  replaced  by  v  ^  , 
fly 

Since  — -  is  zero.  We  obtain: 

8t 

Tr  dVs 

Ps-Vg  .  Vb  —  =  pl'Vs.vl 

ds  ds 


^  -  K3U/ 


(51) 


Up  is  the  relative  velocity  and  is  replaced  by(vs  -  vl).  Furthermore, 
the  expression  *0  .  pL  is  substituted  for  k3,  where  Cr  is  a 

resistance  coefficient  and  is  a  function  of  the  Reynolds  number,  of  the  ' 
particle  with  respect  to  the  surrounding  fluid  and  of  the  acceleration 
parameter.  Since  the  particles  are  always  submerged,  the  Froude 
number  does  not  influence  the  resistance  coefficient.  Introducing  these 


substitutions  and  writing  equation  (51)  fora  unit -volume  we  obtain: 


Pg  .  vs  -  PL  '"L  —  -  Cr  —  PL  (Vs-  vl)  ^ 


dvi 


ds 


ds 


(52) 


The  differential  equation  was  solved  numerically  by  a  step  method .  The 
steps  used  are  ds  =  0.025  m.  At  the  Venturi  entrance  the  solid- 
phase  velocity  vs  is  assumed  equal  to  the  liquid-phase  velocity  vl*  The 
liquid  velocity  vl  can  be  calculated  for  each  cross  section  in  the  Wnturi. 
Since  the  value  Cr  is  not  known  it  has  to  be  assumed.  .This  assumed 
Cr  value  will  lead  to  the  first  approximation.  The  step  method  will  give 
local  values  of  vg  and  so  the  local  relative  velocity  (vg  -  vl)  can  be  com¬ 
puted  .  Furthermore ,  the  average  of  the  relative  velocity  (vg  -  vl)  is 
obtained  and  the  average  change  of  the  relative  velocity  A  (vg  -  vl)  .  We 
shall  now  compare  the  assumed  resistance  coefficient  Cr  with  the  ex¬ 
perimental  results .  Equation  (46)  correlated  the  drag  with  the  loss  due 
to  concentration  and  reads 
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(b-bo)  .  ass  (46) 

A  •  y 

Combining  knoMvn  equation  (46)  and  equation  (41)  we  get; 

1.95  -iv  ‘  c  =  "  (53) 

A*  y 

For  the  drag  equation  we  use  the  same  equation  as  used  in  the  differen¬ 
tial  equation  and  obtain 

1 .95  •  £v  =  .  Cr  (54) 

2g 

Introducing  (vs  -  vl)  as  calculated  with. an  assumed  Cr  value  in  the  first 
approximation,  we. are  able  to  determine  the  assumption  was  satisfactory. 
If  the  assumed  Cr  value  is  satisfactory,  it  wfll  be  taken  for  further  con¬ 
siderations,  if  it  was  not  we  have  to  keep  on  trying,  till  we  find  a  suffi¬ 
ciently  satisfactory  Cr  -  value . 

This  procedure  of  solving  equation  (52)  and  comparing  it  with 
experimental  values  was  done,  due  to  its  elaborate  procedure  for  the 
two  extreme  velocities  in  our  investigation  only .  As  a  minimum  velocity 
2m /s  was  taken  and  4  m/s  was  the  maximum  velocity.  In  Fig.  V-6  the 
velocity  histograph  is  plotted  for  the  extreme  velocities  and  for  both  the 
solid  and  liquid  velocity .  It  was  found  that  the  solid  phase  has  a  smaller 
velocity  at  the  throat  of  the  venturi  and  will  reach  its  maximum  after 
passing  the  throat.  This  can  easily  be  seen  in  the  case  of  the. 2  m/sec 
velocity,  but  is  not  soappmrent  in  the  case  of  4  m/sec,  since  the  chosen 
ds  was  probably  too  big  to  show  this  effect.  But  with  both  velocities  it 
can  be  seen  that,  when  the  particles  leave  the  Venturi  meter  their 
velocity  is  larger  than  that  of  the  liquid . 
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FIGURE  V-6  -VELOCITY  HISTOGRAPH  FOR  LIQUID  AND  SOLID  PHASE 
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The  resulta  of  the  first  and  last  approximation  for  both  extreme 
velocities  are  summarized  in  Table  V-1 . 

In  Fig.  V-7  the  data  are  plotted  on  a  graph  of  Reynolds  number 
versus  drag  coefficient.  In  the  same  graph  the  steady  state  drag  co¬ 
efficients  of  a  sphere,  a  cylinder,  and  a  disk  are  plotted  (Ref.  12).  By 
comparing  our  data  with  the  steady  state  drag  coefficients  ,  we  notice 
that  our  di  ag  coefficients  are  consistently  higher .  Also  on  the  same 
graph  are  shown  drag  coefficients  obtained  with  objects  being  accelerated, 
which  indicate  the  same  trend.  The  research  carried  out  by  Bugliarello 
(Ref.  5)  on  spheres  showed  a  very  similar  effect,  i.e.  at  his  low 
Reynolds  numbers  he  obtained  resistance  coefficients  beyong  Cr  =  10. 
Torobin  &  Gauvin  (Ref.  6)  who  computed  the  acceleration  parameter  for 
Bugliarello 's  results,  ..found  a  variation  of  this  parameter  from  0.01 
to  0.4  .  Bugliarello's  results  are  shown  on. Fig.  V-7.  Laird  et  al. 

(Ref.  7)  obtained  their  data  with  a  decelerating  cylinder  and  observed 
also  this  effect,  that  the  drag  coefficient  was  enlarged.  Keim  (Ref.  8) 
did  his  work  with  an  accelerated  cylinder;  his  data  where  examined  for 
the  acceleration  parameters  by  Torobin  &  Gauvin.  (Ref.  6)  and  are 
shown  in  Fig.,  V-7.  The  data  indicated,  that  an  increase  in  the  .accel¬ 
eration  parameter  is  proportional  to  an  increase  in  the  resistance  co¬ 
efficient.  The  same  phenomenon  can  be  seen  in  the  present  results . 

The  phenomenon  just  described  was  also  noticed  by  Lumon  (Ref.;  .3  and  , 
Ref.  4). 

The  research  of  all  the  references  was  apparently  carried  out 
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in  free  stream  velocity.  In  the  present  work  a  free  stream  vorticity 
exists .  A  paper  dealing  with  the  "Effects  of  fluid  turbulence  on  particle 
drag-coefficient"  by  Torobin  andGsuvin  (Ref.  9)  reports  a  shift  of  the 
laminar -turbulent  boundary  layer  transition  due  to  free  stream  velocity. 
Yet  Ref.  9  does  not  show  in  his  graph  an  observed  deviation  at  other 
Reynolds  numbers  than  the  one  where  the  transition  takes  place.  De¬ 
pending  on  the  intensity  of  the  turbulence  the  transition  is  between 

o  A 

Reynolds  number  and  10  .  Furthermore,  all  of  the  investigations 
conducted  by  the  various  investigators  were  with  smooth  objects .  Though 
it  can  be  said  that  the  roughness  will  be  submerged,  at  these  Reynolds 
numbers ,  in  the  present  tests  in  the  laminar  layer  in  the  case  of  steady 
state,  there  exists  no  information  if  this  is  also  true  in  an  unsteady  case. 


2.2.3  Acceleration  Drag .  We  have  already  mentioned  that  the 
drag  equation  used  in  part  2.2.2  of  this  chapter  should  be  used  where  the 
acceleration  effects  are  of  minor  importance .  With  Cr  10, times  that  of 
uniform  motion  the  acceleration  term  appears  to  be  predominant  and 
equation  (49)  for  the  drag  is  a  poor  choice .  Equation  (48)  would  be  much 
better  .  In  equation  (48)  Up  is  replaced  by(vg  -  vjJ  and  we  have; 

t 

Drag2  =  Ki  (vg  -  +  k  ♦  Vs'Pl^Vs  -  v'l)  (48) 

In  Fig.  V-8  it  is  shown  that  at  each  cross  section  in  the  Venturi  meter, 

(vs  "  vl)2 

the  value  of  the  ratio  - =— -  is  about  the  same.  This  means  that  an 

(vs  -  Vl) 

increase  of  the  second  term  in  equation  (48)  will  cause  a  proportional 


increase  in  the  first  term  .  We  shall  now  write  equation  (48)  for  an  aver¬ 
age  drag  over  the  length  of  the  Venturi  meter.  Then  we  get; 
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Drag2  =  Kj  (V3  -  VL)^  +  k  (Vg-pi)  (vs  -  v'l)  (55) 

We  shall  combine  equation  (55)  with  equation  (53)  and  have 

(1.95-  l^)  .p  =  (Ki  +  k  (Vg  -Pl)  ) ,  c  (56) 

A  .  Y  A*  Y 

In  equation  (56)  the  coefficient  Kj  is  replaced  by  ~  divided 

by  the  concentration  c .  Writing  now  this  equation  for  unit  volume  we  get: 

_  2 

1 .95  .  Kv  =  Cdb  ^  ^  (57) 

2g  g 

(Vg  “  vl)  and  A  (vs  -  vl)  have  been  determined  in  part  2.2.2  of  this 

lyL  V-l 

chapter  and  are  tabulated  in  Fig.  V-8.  From  the  same  table  the 

Reynolds  number  can  be  obtained  and  with  a  given  Reynolds  number  the 

drag  coefficient  Cj)r  can  be  obtained  from  graph  Pig.  V-7.  The  term 

on  the  left  side  of  equation  (57)  is  an  experimental  constant,  already 

discussed  and  g  is  the  acceleration  due  to  gravity .  Everything  is  known  in 

equation  (57) and  it  can  be  solved  for  k  ,  the  added  mass  coefficient. 

The  added  mass  coefficient  was  computed  for  both  extreme  velocities 

and  the  two  different  grain  sizes  and  is  tabulated  in  table  V-l .  It  was 

mentioned  in  part  2.2.1  of  this  chapter,  that  the  added  mass  coefficient 

could  be  corr*3lated  with  a  dimensionless  acceleration  parameter. 

Iversen  and  Balent  (Ref.  2)  were  able  to  find  a  satisfactory  correlation 

for  their  data,  especially  at  high  acceleration  parameters.  Lunnon 

(Ref.  3  and  Ref.  4)  reported  to  have  obtained  values  of  k  ranging  from 

0.5  to  2.0.  Bugliarello  (Ref.  5)  failed  to  give  a  satisfactory  correlation, 

neither 'did  his  .values  correlate  with  ,  nor  were  they  in  the  range 

V2 

reported  by  Lunnon.  The  two  last  authors  investigated  a  sphere,  while 


FIGURE  V-9  -  VIRTUAL  MASS  COEFFICIENT 


Iveraen  and  Balent  did  their  work  on  a  disk.  In  Fig.  V-9  the  authors' 
data  are  compared  with  the  references  just  discussed .  One  might  see 
that  our  data  fall  at  least  in  a  reasormble  range  (0.45  -  k  -  1 .70  ) . 

Since  our  scope  was  not  to  investigate  the  validity  of  the  added  mass 
coefficient,  nor  its  dependence  from  the  acceleration  parameter,  the 
number  of  values  obtained  is  limited.  The  particle  grain  size  range 
and  the  velocity  range  were  just  too  small .  All  we  can  conclude  is  that 
the  magnitude  of  the  added  mass  coefficient  is  reasonable  and  provided 
probably  a  possible  explanation  of  the  high  drag  or  the  large  energy 


loss . 


V 


Table  V-1 
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VI.  SUMMARY.  CONCLUSIONS,  AND  SUGGESTIONS 

The  investigation  described  in  this  report  was  concerned  with  two- 
phase  flow  in  a  closed  pipe  system. 

A  device  which  gives  information  on  flow  rate  of  mixture  and  the 
concentration  was  designed  ,  tested,  and  found  useful.  This  device, 
called  a  "Loop",  consists  of  two  vertical  pipes,  with  opposite  flow  di¬ 
rections.  Combination. (the  summation  and  the  difference)  of  the  two 
head  readings  for  each  flow  gave  us  the  desired  resiilts.  The  effect  was 
explained  by  a  theory . 

Fxirthermore ,  a  Venturi  meter  was  tested  for  its  applicability  in 
a  two-phase  flow.  Besides  the  usually  recorded  pressure  drop,  the 
pressure  loss  across  the  horizontal  Venturi  meter  was  obtained.  It  is 
shown  by  experiment  and  theory,  that  the  pressime  readings  in  terms  of 
mixture  gives  information  of  the  flow  rate.  The  energy  loss  was  corre¬ 
lated  with  the  concentration .  An  attempt  was  ma  de  to  explain  also  this 
by  a  theory.  Since  the  concept  of  virtiial  mats  seems  to  be  important, 
and  the  existing  literature  is  quite  incomplete,  our  explanation  is  of 
limited  value . 

For  all  practical  purposes  it  can  be  stated  that  the  "Loop" 
system  will  provide. a  possibility  to  determine  immediately  the  flow  rate 
of  mixture  and  the  concentration,  if  two  manometer  readings  one  at  a  riser 
section  and  one  on  downcomer  section  are  known.  The  practical  appli¬ 
cation  of  the  modified  Venturi  meter  is  not  quite  as  straight  forward. 


Il- 
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In  order  to  determine  the  concentration. the  flow  rate  must  be  known,  but 
in  order  to  obtain  the  flow  rate,  the  concentration  has  to  be  known.  But 
fast  trial  and  error  solution  will  lead  to  a  reliable  determination  of  flow 
rate  and  concentration. 

The  two-phase  flow  problem  reqtiires  additional  research.  In  a 
rather  extensive  way  straight  pipe  tests  were  already  investigated.  But 
not  too  much  has  been  done  in  studying  different  units  in  pipeline  systems, 
their  applicabilities  and  their  effects.  The  present  study  was  concerned 

with  a  flow  meter,  the  Ventiuri  tube  in  a  horizontal  position.  An.immedi- 

1 

ate  suggestion  could  be  to  investigate  this  same  flow  meter  in  vertical  or 
oblique  pipe  systems Also  it  would  be  of  interest  to  study  the  behavior 
of  other  flow  meters  like  the  orifice  meter,  metering  nozzle  and  elbow 
meters ,  in  a  sand-water  mixture . 

It  is  suggested  therefore  that  a  study  with  sand  of  different  ge¬ 
ometry  and  density  would  make  it  possible  to  establish  more  general 
constants .  It  is  believed  that  the  pipe  diameter  will  not  change  the  re¬ 
sults,  although  no  attempt  was  made  to  check  this. 

The  results  obtained  at  the  Venturi  meter  can  probably  be  plotted 
in  form  of  a  nomograph ,  which  would  provide  a  fast  interpretation  of  the 
Venturi  meter  reading. 
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SYMBOL 

A 

Al 

A2 

am 

AD 

W 

b 

b  -  bo 

CD 

CR 


(1-c) 

Cdr 

Cr 

Cy 

d50 

da 

Dragi 

Pragg 


VIII.  NOMENCLATURE 

EXPLANATION 

Area 

Area  of  pipe  cross  section 

Area  of  cross  sevjtion  at  throat  of  Venturi-meter 

Pressure  drop  in  terms  of  mixture 

Acceleration  modulus:  ^celeration  x  Length 

(Velocity)^ 

total  energy  loss  across  Venturi  meter 

energy  loss. across  Venturi  meter,  due  to  concentration 

concentration  in  downcomer  section 

concentration  in  riser  section 

average  concentration  in  volume  percentage;  taken  up  by 
splid  .phase 

average  concentration  in  volume  percentage;  taken  up 

by  liquid  phase . 

steady  state  drag  coefficient 

resistance  coefficient ,  dependent  on  steady  state  and 
and  unsteady  state  effects 

coefficient  in  flow  rate  equation  of  Venturi  meter,  taking 

care  of  real  fluid  effects 

average  grain  size  diameter 

or  As  stepwise  increment  in  differ 

steady  state  drag 

summation  of  steady  state  and  unsteady  state  drag 


SYMBOL 


EXPLANATION 
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Drags 

F 

Ffr 

Fp 

g 

L 

iy 

k 

Kl 

K3 

m 

Q 

Qs 

Ql 

Rd 

Rr 

Rfr 

Rsuap 

Up 

Vs 

V  or  vm 

Vs 

VL 


combined  steady  state  and  unsteady  state  drag 
force 

wall  friction. force 
particle  force 
acceledation  due  to  gravity 
length  of  riser  and  downcomer  section 
length  over  which  energy  loss  at.  Venturi  meter  is  re¬ 
corded 

added  or  virtual  mass  coefficient 
steady  state  constant,  including  Cdj^ 
steady  and  unsteady  state  constant 
symbol  for  mass 

or  Qm  flow  rate  of  mixture;  Q  =  Qs  +  Ql 
flow  rate  of  solid  phase 
flow  rate  of  liquid  phase 
manometer  reading  at  downcomer  section 
manometer  reading  at  riser  section 
manometer  reading  due  to  friction 
manometer  reading  due  to  suspension 
relative  velocity  of  liquid  and  solid  phase  (Vp  ) 
volume  of  solid  particle 
or  velocity  of  mixture,  m  =  Q/A 
velocity  of  solid  phase 
velocity  of  liquid  phase 
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SYMBOL 

EXPLANATION 

^settl 

settling  velocity 

a 

ratio  of  flow  rate  of  solids  to  flow  rate  of  mixture 

1  -  O' 

ratio  of  flow  rate  of  liquids  to  flow  rate  of  mixture 

yl 

specific  weight  of  liquid  phase 

Ts 

specific  weight  of  solid  phase 

Ym 

specific  weight  of  mixture 

Pl 

density  of  liquid  phase 

Pb 

density  of  solid  phase 

Pm 

density  of  mixtiare 
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Ref.  1 

.  Ref.  2 

Ref.  3 
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Ref.  6 
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Summary  of  Data  Loop  System 


1 

2 

3 

4 

5 

6 

7 

8 

flow  rate 

sum. 

V  niixt. 

sum. 

Manometer  Readings 
riser  downcomer 

diff . 

cone . 

no. 

L/sec 

m 

m 

m 

m 

m 

vol.  % 

series 

. 

run 

- 

- 

4+5 

-  ■ 

- 

4-5 

- 

Clear  Water 

10.00* 

0.142 

0.142 

0.071 

0.071 

0 

0 

19 

10.55* 

0.152 

0.152 

0.076 

0.076 

0 

0 

20 

12.10* 

0.214 

0.214 

0.107 

0.107 

0 

0 

30 

13.82* 

0.292 

0.292 

0.146 

0.146 

0 

0 

40 

15.00* 

0.348 

0.348 

0.175 

0.173 

0 

0 

50 

16.55* 

0.388 

0.388 

0.194 

0.194 

0 

0 

60 

*  An. average  over 

a  number  of  measurements  is  recorded 

0 

7.7 

0.113 

0.143 

Sand  No.  2 

0.338  -0.195 

0..533 

16.5 

19-1 

8.6 

0.125 

0.140 

0.323 

-0.183 

0.406 

8* 

2 

8.4 

0.122 

0.144 

0.318 

-0.174 

0.492 

11.4 

3 

9.8 

0.137 

0.140 

0.305 

-0.165 

0.470 

3.8 

4 

8.8 

0.116 

0.138 

0.290 

-0.152 

0.442 

11.0 

5 

9.1 

0.125 

0.143 

0.268 

-0.125 

0.393 

8.5 

6 

9.1 

0.125 

0.144 

0.244 

-0.100 

0.344 

8.5 

7 

9.7 

0.131 

0.147 

0.223 

-0.076 

0.299 

6.1 

8 

9.4 

0.131 

0.143 

0.201 

-0.058 

0.259 

5.3 

9 

8.9 

0.128 

0.143 

0.177 

-0.034 

0.201 

6.2 

10 

9.35 

0.131 

0.140 

0.155 

-0.015 

0.170 

4.5 

11 

9.70 

0.125 

0.131 

0.137 

-0.006 

0.143 

3.4 

12 

10.5 

0.125 

0.128 

0.122 

+0.006 

0.116 

1..3 

13 

9.4 

0.125 

0.131 

0.113 

+0.018 

0.095 

2:4 

14 

9.80 

0.128 

0.134 

0.107 

+0.027 

0.080 

2.3 

15 

10.05 

0.131 

0.131 

0.098 

+0.033 

0.065 

0.9 

16 

10.7 

0.131 

0.131 

0.091 

+0.040 

0.051 

0.8 

17 

8.5 

0.119 

0.140 

0.280 

-0.140 

0.420 

10.8 

19-18 

9.3 

0.122 

0.137 

0.259 

-0.122 

0.381 

7.4 

19 

9.7 

0.125 

0.140 

0.241 

-0.101 

0.342 

7.0 

20 

9.2 

0.128 

0.144 

0.223 

-0.079 

0.302 

8.4 

21 

9.4 

0.131 

0.143 

0.207 

-0.064 

0.271 

5.4 

22 

9.3 

0.125 

0.131 

0.183 

-0.052 

0.235 

3.5 

23 

9.9 

0.122 

0.132 

0.162 

-0.030 

0.192 

4.2 

24 

9.9 

0.125 

0.135 

0.147 

-0.012 

0.159 

3.7 

25 
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1 

flow  rate 

2  3 

sum. 

y  mixt .  sum  . 

4  5  6 

Manometer  Eeadings 
riser  downcomer 

7 

cone . 

8 

no 

• 

L/sec 

m 

m 

m 

m 

m 

vol.  % 

series 

run 

- 

- 

4+5 

- 

- 

4-5 

- 

_ 

Sand  No .  2 

9.3 

0.122 

0.131 

0.134 

-0.003 

0.137 

4.4 

19- 

■26 

9.5 

0.128 

0.134 

0.125 

+0.009 

0.116 

1.9 

27 

9.5 

0.128 

0.134 

0.116 

+0.018 

0.098 

3.3 

28 

10.45 

0.128 

0.134 

0.167 

+0.027 

0.080 

1.0 

29 

8.45 

0.119 

0.150 

0.304 

-0.154 

0.458 

11.4 

19- 

-30 

9.00 

0.119 

0.148 

0.286 

-0.138 

0.424 

9.6 

31 

9.35 

0.125 

0.145 

0.263 

-0.118 

0.381 

6.4 

32 

8.60 

0.119 

0.141 

0.244 

-0.103 

0.347 

9.0 

33 

9.60 

0.125 

0.140 

0.226 

-0.086 

0.302 

4.8 

34 

9.10 

0.128 

0.146 

0.207 

-0.061 

0.268 

6.0 

35 

9.4 

0.131 

0.148 

0.194 

-0.046 

0.240 

4.9 

36 

10.00 

,0.131 

0.142 

0.172 

-0.030 

0.202 

2.7 

37 

9.7 

0.131 

0.147 

0.157 

-0.015 

0.172 

3.3 

38 

9.9 

0.134 

0.140 

0.147 

-0.003 

0.150 

3.3 

39 

10.45 

0.134 

0.145 

0.135 

+0.010 

0.125 

3.1 

40 

9.9 

0.137 

0.146 

0.125 

+0.021 

0.104 

2.8 

41 

9.8 

0.137 

0.143 

0.113 

+0.030 

0.083 

1.4 

42 

10.45 

0.137 

0.143 

0.104 

+0.039 

0.065 

1.3 

43 

10.40 

0.137 

0.143 

0.097 

+0.046 

0.051 

0.9 

44 

10,00 

0.134 

0.137 

0.091 

+0.046 

0.045 

0.5 

45 

9.30 

0.134 

0.136 

0.088 

+0.048 

0.040 

0.9 

46 

7.0 

0.098 

0.132 

0.360 

-0.228 

0.588 

17.6 

19 

.47 

8.05 

0.107 

0.139 

0.358 

-0.219 

0.577 

14.2 

48 

8.4 

0.122 

0.148 

0.352 

-0.204 

0.556 

9.2 

49 

8.4 

0.119 

0.150 

0.339 

-0.189 

0.528 

10.4 

50 

8.7 

0.125 

0.156 

0.326 

-0.170 

0.496 

10.5 

51 

9.5 

0.131 

0.156 

0.308 

-0.152 

0.452 

7.7 

52 

9.3 

0.125 

0.151 

0.282 

-0.131 

0.413 

9.3 

53 

8.8 

0.125 

0.150 

0.260 

-0.110 

0.370 

9.3 

54 

9.4 

0.125 

0.145 

0.242 

-0.097 

0.339 

6.2 

55 

9.6 

0.125 

0.144 

0.220 

-0.076 

0.396 

6.2 

56 

9.3 

0.131 

0.146 

0.201 

-0.055 

0.256 

3.9 

57 

10.1 

0.131 

0.142 

0.188 

-0.046 

0.234 

2.7 

58 

10.4 

0.128 

0.143 

0.173 

-0.030 

0.203 

4.0 

59 

10.4 

0.134 

0.145 

0.160 

-0.015 

0.175 

3.1 

60 

10.1 

0.134 

0.145 

0.151 

-0.006 

0.157 

2.2 

61 
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1 

2 

3 

4 

5 

6 

7 

8 

sum . 

Manometer  Readings 

diff. 

cone . 

no. 

flow  rate 

y  mixt. 

sum . 

riser 

downcomer 

L/sec 

m 

m 

m 

m 

m 

vol .% 

series 

run 

- 

- 

4+5 

- 

- 

4-5 

- 

- 

Sand  No .  2 

11.6 

0.140 

0.144 

0.138 

+0.006 

0..13.2 

,0.8 

19-62 

9.7 

0.131 

0.141 

0.126 

+0.015 

0.111 

2.8 

63 

9.8 

0.137 

0.140 

0.116 

+0.024 

0.092 

0.9 

64 

9.7 

0.134 

0.143 

0.113 

+0.030 

0.083 

2.8 

65 

9.9 

0.140 

0.147 

0.107 

+0.040 

0.067 

0.5 

66 

10.6 

0.140 

0.140 

0.094 

+0.046 

0.048 

0.9 

67 

i:o..,o 

0.137 

0.137 

0.088 

+0.049 

0.039 

0.5 

68 

:lo..4 

0.134 

0.137 

0.085 

+0.052 

0.032 

0.5 

69 

8.2 

0.155 

0.204 

0.430 

-0.226 

0.656 

1.8 

30-1 

8.5 

0.155 

0.205 

0,404 

-0.199 

0.603 

1.9 

2 

10.2 

0.189 

0.206 

0.386 

-0.180 

0.566 

6.7 

3 

9.6 

0.163 

0.208 

0.367 

-0.159 

0.526 

17.5 

4 

10.3 

0.180 

0.214 

0.345 

-0.131 

0.476 

11.9 

5 

78 


Summary 

of  Data 
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1 

2 

3 

4 

5 

6 

7 

8 

sum . 

Manometer  Readings 

flow  rate 

y  mixt 

.  sum . 

riser 

downcomer 

diff. 

cone . 

no. 

L/sec 

m 

m 

m 

m 

m 

vol  % 

series 

run 

- 

- 

- 

4+5 

- 

4-5 

- 

- 

Sand  No.  2 

11  .2 

0.189 

0.217 

0.318 

-0.101 

0.419 

8.9 

30-6 

10.15 

0.186 

0.217 

0.294 

-0. 077 

0.371 

10.3 

7 

10.8 

0.165 

0.212 

0.266 

-0.058 

0.324 

4.2 

8 

10.5 

0.189 

0.210 

0.238 

-0.028 

0.266 

7.3 

9 

11.7 

0.204 

0.214 

0.214 

-0.000 

0.214 

2.3 

10 

11  .3 

0.210 

0.220 

0.199 

+0.021 

0.178 

2.4 

11 

11 .7 

0.204 

0.221 

0.183 

+0.038 

0.145 

5.1 

12 

12.0 

0.204 

0.214 

0.162 

+0.052 

0.110 

2.7 

13 

12.6 

0.205 

0.217 

0.150 

+0.067 

0.083 

3.3 

14 

12.4 

0.210 

0.214 

0.138 

+0.076 

0.062 

1.5 

15 

10.4 

0.207 

0.213 

0.134 

+0.079 

0..055 

2.2 

16 

11.5 

0.204 

0.207 

0.128 

+0.079 

0.049 

1.2 

17 

12.5 

0.204 

0.207 

0.122 

+0.085 

0.037 

0.7 

18 

8.5 

0.159 

0.197 

0.425 

-0.228 

0.653 

14.3 

30-19 

9.1 

0.165 

0.202 

0.412 

-0.210 

0.622 

14.5 

20 

9.7 

0.155 

0.204 

0.388 

-0.192 

0.580 

15.5- 

21 

10.7 

0.178 

0.207 

0.368 

-0.162 

0.530 

8.9 

22 

10.5 

0.186 

0.210 

0.330 

-0.120 

0.450 

8.2 

23 

10.5 

0.186 

0.209 

0.300 

-0.091 

0.391 

9.1 

24 

11 .1 

0.186 

0.212 

0.276 

-0.064 

0.340 

8.6 

25 

14.6 

0.201 

0.224 

0.256 

-0.052 

0.288 

5.6 

26 

12.0 

0.216 

0.226 

0.232 

-0.006 

0.238 

2.6 

27 

11 .3 

0.213 

0.203 

0.214 

+0.009 

0.205 

3.2 

28 

12.2 

0.198 

0.210 

0.192 

+0.018 

0.174 

4.1 

29 

11.3 

0.186 

0.196 

0.155 

+0.041 

0.114 

2.8 

30 

12.5 

0.186 

0.193 

0.137 

+0.056 

0.181 

2.6 

31 

12.4 

0.195 

0.205 

0.140 

+0.065 

0.075 

1.1 

32 

12.3 

0.202 

0.206 

0.128 

+0.078 

0.050 

0,7 

33 

12.0 

0.202 

0.205 

0.122 

+0.083 

0.039 

0.8 

34 

8.2 

0.151 

0.208 

0.424 

-0.216 

0.640 

22.0 

30-40 

9.4 

0.171 

0.202 

0.400 

-0.198 

0.598 

12.1 

41 

11.0 

0.164 

0.195 

0.378 

-0.183 

0.561 

11.2 

42 

10.1 

0.161 

0.197 

0v360 

-0.163 

0.523 

13.5 

43 

12.0 

0.183 

0.204 

0.335 

-0.131 

0.466 

6.9 

44 

9.9 

0.164 

0.203 

0,310 

-0.107 

0.417 

15.3 

45 

12.1 

0.183 

0.204 

0.284 

-0.080 

0.364 

7.2 

46 

12.4 

0.189 

0.209 

0.258 

-0.049 

0.307 

6.6 

47 

12.4 

0.208 

0.213 

0.240 

-0.027 

0.267 

1.8 

48 

11.9 

0.192 

0.213 

0.222 

-0.009 

0.231 

7.3 

49 
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1 

2 

3 

4 

5 

6 

7 

8 

flow  rate 

sum . 
y  mixt. 

sum . 

Manometer  Readings 
riser  downcomer 

diff. 

cone . 

no. 

L/sec 

m 

m 

m 

m 

m 

vol.% 

series 

nm 

- 

- 

- 

4+5 

- 

4-5 

- 

- 

Sand  No.  2 

12.3 

0.208 

0.187 

0.202 

+0.005 

0.187 

2.9 

30-50 

7.7 

0.157 

0.203 

0.413 

-0.210 

0.413 

19.5 

30-51 

10.2 

0.171 

0.203 

0.398 

-0.195 

0.398 

12.1 

52 

10.6 

0.174 

0.208 

0.385 

-0.177 

0.385 

12.5 

53 

10.1 

0.177 

0.192 

0.364 

-0.152 

0.344 

12.2 

54 

10.5 

0.168 

0.202 

0.336 

-0.134 

0.336 

7.8 

55 

11.7 

0.183 

0.210 

0.314 

-0.104 

0.314 

8.6 

56 

12.0 

0.204 

0.213 

0.293 

-0.080 

0.293 

5.0 

57 

11.4 

0.189 

0.210 

0.268 

-0.058 

0.268 

6.8 

58 

11.8 

0.192 

0.211 

0.247 

-0.036 

0.247 

5.8 

59 

12.5 

0.198 

0.214 

0.232 

-0.018 

0.232 

4.4 

60 

11.4 

0.189 

0.203 

0.214 

-0.009 

0.212 

5.2 

61 

11.8 

0.180 

0.196 

0.190 

+0.006 

0.184 

5.8 

62 

12.3 

0.192 

0.204 

0.174 

+0.030 

0.144 

3.7 

63 

11.4 

0.195 

0.205 

0.165 

+0.040 

0.120 

3.2 

64 

12.1 

0.198 

0.204 

0.155 

+0.049 

0.106 

2.6 

65 

11.9 

0.198 

0.208 

0.150 

+0.058 

0.092 

2.3 

66 

12.0 

0.201 

0.204 

0.140 

+0.064 

0.076 

0.8 

67 

12.2 

0.198 

0.204 

0.131 

+0.073 

0.058 

1.1 

68 

11.7 

0.198 

0.201 

0.128 

+0.073 

0.055 

0.8 

69 

11 .6 

0.198 

0.201 

0.122 

+0.079 

0.043 

0.8 

70 

11.2 

0.201 

0.24 

0.50 

-0.256 

0.760 

12.5 

40-1 

11.4 

0.204 

0.25 

0.48 

-0.229 

0.715 

13.0 

2 

11.8 

0.210 

0.25 

0.44 

-0.192 

0.63 

13.4 

3 

11.6 

0.216 

0.26 

0.42 

-0.159 

0.58 

14.0 

4 

11.7 

0.229 

0.27 

0.38 

-0.107 

0.49 

13.1 

5 

12.1 

0.253 

0.29 

0.33 

-0.043 

0.37 

7.8 

6 

13.4 

0.259 

0.29 

0.28 

+0.006 

0.29 

7.4 

7 

13.6 

0.265 

0.28 

0.23 

+0.05 

0.18 

4.9 

8 

13.8 

0.277 

0.28 

0.20 

+0.08 

0.12 

1.3 

9 

14.1 

0.277 

0.27 

0.18 

+0.09 

0.09 

0.9 

10 

13.9 

0.277 

0.27 

0.17 

+0.107 

0.06 

1.0 

11 

10.1 

0.201 

0.25 

0.48 

-0.23 

0.71 

17.5 

40-12 

11.2 

0.208 

0.26 

0.45 

-0.19 

0.64 

15.8 

13 

- 

- 

- 

0.42 

-0.15 

0.57 

15.1 

14 

12.3 

0.238 

0.28 

0.39 

-0.11 

0.50 

11.4 

15 

13.4 

0.275 

0.26 

0.31 

0.00 

0.31 

7.1 

17 

12.4 

0.289 

0.29 

0.26 

+0.03 

0.23 

3.0 

18 

ft' 
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3 


Loop  System 
6  7  8 


80 


1 


4 


sum .  Manometer  Reading 


flow  rate 

y  mixt 

.  sum. 

riser 

downcomer 

diff. 

cone . 

no. 

L/sec 

m 

m 

m 

m 

m 

vol.% 

series 

run 

- 

- 

- 

4+5 

- 

4-5 

- 

- 

Sand  No .  2 

13.4 

0.259 

0.28 

0.21 

+0.07 

0.14 

4.10 

40-19 

13.7 

0.275 

0.28 

0.18 

+0.10 

0.08 

1.30 

20 

13.8 

0.268 

0.28 

0.17 

+0.11 

0.06 

1.70 

21 

14.2 

0.275 

0.27 

0.16 

+0.11 

0.04 

0.60 

22 

14.0 

0.278 

0.27 

0.15 

+0.12 

0.03 

0.96 

23 

14.0 

0.278 

0.28 

0.15 

+0.13 

0.02 

0.30 

24 

10.2 

0.217 

0.28 

0.42 

-0.14 

0.56 

16.0 

40-25 

12.0 

0...228 

0.28 

0.39 

-0.11 

0.50 

12.5 

26 

12. 5 

0.240 

0.27 

0.34 

-0.07 

0.41 

9.6 

27 

11.7 

0.253 

0.29 

0.31 

-0.02 

0.33 

6.1 

28 

12.7 

0.292 

0.31 

0.29 

+0.02 

0.27 

2.5 

29 

10.4 

0.210 

0.28 

0 .45 

-0.17 

0.62 

17.6 

40-30 

13.0 

0.222 

0.27 

0.43 

-0.16 

0.59 

12.3 

31 

12.3 

0.226 

0.28 

0.40 

-0.12 

0.52 

13.0 

32 

12.5 

0.243 

0.29 

0.36 

-0.07 

0.43 

9.5 

33 

12.9 

0.258 

0.28 

0.32 

-0.04 

0.36 

5.7 

34 

12.4 

0.254 

0.29 

0.29 

-0.00 

0.29 

8.8 

35 

14.6 

0.268 

0.29 

0.25 

+0.04 

0.21 

5.0 

36 

11.3 

0.222 

0.28 

0.42 

-0.14 

0.54 

13.7 

40-37 

12.7 

0.228 

0.27 

0.38 

-0.11 

0.49 

10.4 

.  38 

12.6 

0.250 

0.25 

0.30 

-0.05 

0.35 

6.8 

40 

13.5 

0.256 

0.25 

0.27 

+0.02 

0.25 

7.5 

41 

13.6 

0.265 

0.29 

0.24 

+0.05 

0.19 

5.4 

42 

13.1 

0.271 

0.28 

0.22 

+0.06 

0.16 

1.05 

43 

11.8 

0.226 

0.26 

0.37 

-0.11 

0.48 

12.0 

40-44 

13.0 

0.235 

0.27 

0.32 

-0.05 

0.37 

8.8 

45 

13.4 

0.253 

0.27 

0.29 

-0.02 

0.31 

5.1 

46 

13.1 

0.253 

0.27 

0.27 

+0.00 

0.27 

3.8 

47 

12.5 

0.259 

0.28 

0.26 

+0.02 

0.24 

4.4 

48 

12.5 

0.268 

0.28 

0.24 

+0.04 

0.20 

2.8 

49 

13.8 

0.256 

0.28 
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